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1.0  INTRODUCTION 


As  radar  systems  designs  become  more  complex  and  versatile,  their  per¬ 
formance  becomes  Increasingly  sensitive  to  the  operational  environment.  This 
places  an  increased  burden  on  the  designer  to  incorporate  the  effects  of  the 
environment  in  system  design  studies.  However,  before  this  can  be  done,  it  is 
necessary  to  develop  an  accurate  model  of  the  environment.  In  the  case  of 
ground  clutter  or  multipath,  this  means  that  there  is  a  need  for  rough  sur- 
face  scattering  models.  Such  models  must  not  only  be  based  upon  sound  physi¬ 
cal  principles  but  also  exhibit  agreement  with  measurements. 

The  purpose  of  this  study  is  to  provide  improved  models  for  surface  scat¬ 
tering  from  terrain  in  the  microwave  frequency  range  and  near  grazing  inci¬ 
dence.  The  basic  approach  entails  applying  the  composite  surface  scattering 
theory  to  a  lossy,  randomly  rough,  dielectric  surface.  As  long  as  the  inci¬ 
dent  energy  is  not  too  near  grazing  and  the  surface  is  reasonably  free  of 
sharp  edges  or  cusps,  the  composite  model  should  be  a  reasonable  description 
of  the  scattering  process.  Since  the  composite  model  is  based  upon  the  combin¬ 
ing  of  two  asymptotic  scattering  theories,  it  is  approximate.  Thus,  an  addi¬ 
tional  goal  of  this  study  is  to  investigate  new  techniques  for  improving  the 
composite  model.  As  our  understanding  of  the  scattering  process  increases 
and  the  models  become  more  accurate,  complex  factors  such  as  vegetation  and 
snow  cover  can  be  studied,  analyzed,  and  incorporated  into  the  basic  model. 
Thus,  the  overall  effort  comprises  a  phased  approach. 

This  interim  report  details  the  initial  results  of  the  study.  Present 
and  future  efforts  on  this  investigation  will  be  concerned  with  a  more  detailed 
study  of  some  of  the  results  reported  here  and  also  obtaining  improved  models 
for  the  coherent  or  average  scattered  field. 


1.1  Sunnary  of  Results 


Section  2  corrects  an  error  in  the  basic  composite  scattering  model. 

In  particular,  it  is  shown  that  shadowing  is  improperly  accounted  for  and  this 
error  is  corrected.  In  the  corrected  version  of  the  composite  mqdel,  it  is 
shown  that  the  conventional  shadowing  function  multiplies  both  the  zeroth  and 
first  order  incoherent  scattered  power  perturbation  terms.  Thus,  even  for  a 
perfectly  conducting  surface  the  first  order  term  will  go  to  zero  for  back- 
scattering  at  grazing  incidence  due  to  the  shadowing  function.  The  first 
order  perturbation  power  suffers  an  additional  attenuation  due  to  the  shadow¬ 
ing  of  unfavorably  oriented  large  scale  surface  slopes;  however,  this  effect 
is  relatively  small  compared  to  the  impact  of  the  conventional  shadowing  func¬ 
tion. 

Section  3  demonstrates  how  the  shadowing  function  for  non-Gaussian  may 
be  quite  easily  obtained  from  existing  shadowing  theories.  The  important 
surface  characteristic  in  the  general  case  is  the  probability  density  func¬ 
tion  of  the  large  scale  slopes  in  the  plane  of  incidence.  Exp licit  results 
are  obtained  for  a  surface  characterized  by  roughness  whose  probability  den¬ 
sity  function  is  exponential.  The  results  of  this  study  are  particularly  im¬ 
portant  for  terrain  scattering  because  terrain  cannot  always  be  described  by 
a  Gaussian  height  probability  density  function. 

Section  4  extends  the  composite  model  to  bistatic  scattering  from  a 
lossy,  dielectric,  rough  surface.  The  details  are  presented  for  three  cases 
of  increasing  complexity;  backscattering  from  a  surface  with  only  small  scale 
roughness,  biatatlc  scattering  from  a  surface  with  only  small  scale  roughness, 
and,  finally,  bistatic  scattering  from  a  composite  (large  and  small  scales  of 
roughness)  surface.  This  approach  is  a  logical  progression  from  the  simple 
to  the  complex  and  is  therefore  beneficial  to  the  reader.  Furthermore,  this 
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approach  facilitates  checking  the  results  of  the  perturbation  theory  against 
existing  solutions. 

Section  5  discusses  one  technique  for  relating  the  joint  probability 
density  function  of  the  surface  heights  to  the  joint  density  function  for  the 
slopes.  The  technique  was  originally  obtained  from  an  analysis  of  optical 
scattering  from  a  rough  surface  but  its  relevance  to  this  problem  has  apparently 
been  overlooked.  Although  the  technique  is  not  always  applicable  to  measured 
data,  there  are  cases  where  it  can  provide  the  desired  transformation. 
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2.0  A  CORRECTION  TO  THE  COMPOSITE  SURFACE  SCATTERING  MODEL 

2 . 1  Background 

In  [1],  a  solution  to  the  problem  of  backscattering  from  a  randomly  rough, 
perfectly  conducting  surface  comprising  both  large  and  small  scales  of  rough¬ 
ness  was  presented.  As  a  direct  consequence  of  the  stipulation  that  the  sur¬ 
face  roughness  was  a  zero  mean  jointly  Gaussian  process,  the  scattering  cross 

section  per  unit  area  was  determined  to  be  the  sum  of  two  terms,  i.e. 

*.  _ . 

app.(9,<J>)  -  fap°t(8,<{,)3o  +  £crpp »  C©»4>) 3-l  •  The  [o  °» (8,<j>)]0  contribution  is 
dominant  near  normal  incidence  and  results  from  the  shadow  corrected  optical 
like  reflection  from  properly  oriented  facets  or  specular  points  on  the  sur¬ 
face.  The  [Op°,  (9  ,(}>)  term  is  due  to  Bragg  resonance  scattering  from  the 
small  scale  surface  features  with  appropriate  accounting  for  the  resonance 
broadening  effects  of  the  large  scale  surface  undulations.  Although  shadowing 
is  formally  accounted  for  in  a  correct  manner  in  the  [Op°,  (8  ,<|>)  ]^  term, 
there  is  an  error  In  the  exact  representation  of  the  shadowing  function  which 
leads  to  an  incorrect  estimate  of  the  effects  of  shadowing  on  large  angle  of 
incidence  scattering.  The  goal  of  this  section  is  to  correct  the  above  error 
and  to  properly  account  for  the  effect  of  large  scale  shadowing  on  the  small 
scale  scattering  term. 

2.2  Discussion  of  the  Error 

The  analysis  presented  in  [1]  relating  to  the  determination  of  [opp,  (8,<J>)  ^ 
is  correct  up  to  and  including  equation  (23).  The  problem  with  the  analysis 
following  (23)  is  a  result  of  inadequate  attention  to  the  definition  of  the 
shadowing  function  R(8,4>)  .  That  is,  R(8,<J>)  ,  as  it  appears  in  (23)  of  [1], 
can  and  should  be  expressed  as  follows; 

M8.W  -  |  Cju  -  V  ,  %  -  V)  G-” 
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where 

P<  *(*1  |  hx  ’hy)  “  <p(  |  Zl  •  Qlx  *  ^Ay)/^  (2-2) 

is  the  probability  that  an  Incident  ray  having  direction  will  intersect  a 
point  on  the  large  scale  surface  with  orthogonal  slopes  £ and  and  will 
not  be  shadowed  by  any  other  part  of  the  surface  regardless  of  the  height  £0 


of  the  point  in  question.  The  symbol  <•>  in  (2.2)  denotes  the  ensemble 

average  over  all  values  of  large  scale  height  £«  and  the  notation  P^(») 

is  the  same  as  that  employed  by  Sancer  [2 ]  in  his  excellent  analysis  of  the 

effect  of  shadowing  on  [Op°i  (9»<J>)]g  •  It  should  be  noted  that  as  a  result 

of  the  integrations  in  (23)  of  [1],  the  shadowing  function  R(0,<f>)  is  equal 

to  P^°\*)  evaluated  at  the  specific  large  scale  slope  values  given  by 

C0  -  k  /B  and  CD  ■  k  /B  where  k  *  -2k  sin  0  cos  4>  ,  k  *  -2k  sin  0  sin  4>  , 
lx  ox  ly  oy  ox  o  ’  oy  o 

B  *  2kQcos  6,9  is  the  angle  of  incidence  relative  to  the  normal  to  the  mean 
(C  -0)  plane,  and  <J>  is  the  azimuth  direction  of  incidence. 

In  order  to  more  clearly  understand  and  therefore  rectify  the  error  in 
[1],  it  is  beneficial  to  repeat  equation  (23)  of  [1],  i.e. 


00  00 

J  J  ^I(xl,yl)I(x2,y2)fpp,(^£x1,CAx2’CAy1,CAy2)exp[JBCCA2"CA1)]^ 


— OO  —CO 


/  K  k  k  ^  \ 

eXp(_JkoxAx"JkoyAy)  dAxdAy  3  R(e’<{,)fpp,\‘^,~B1,"BJL,"BiL/ 


_1_ 

2tt 


j *  J  exp  |  -  j  (kQxAx  +  kQyAy)  -  4k^  cos2  0  [l  -  P^(Ax,Ay)]  j  dAxdAy 


(2.3) 
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where  Ax  ■  -  x^  and  Ay^y^-Yi  .  In  (2.3)  iAc^  ,  y.^  is  one  if  the  point 

^2  2 

(X1  *  yl  *  )  on  the  -*-arSe  scale  surface  having  slopes  and  ?£y  is 

2  2  2  2  2 

illuminated  and  zero  if  it  is  shadowed.  The  function  f  ,  is  defined  in  [1], 

PP 

2 

is  the  mean  square  height  of  the  large  scale  surface,  and  p^(*)  is  the 

normalized  autocorrelation  function  of  the  large  scale  surface  height.  It 

should  be  noted  that  the  left  hand  side  of  (2.3)  represents  the  Fourier  trans- 

form  of  the  <•>  term  (from  AxA  y -  space  to  k  k  -space)  and  (2.3)  is 

r  ox  oy 

supposed  to  show  where  the  transform  variables  appear  in  the  result.  Such  know¬ 
ledge  is  essential  to  accomplishing  the  convolution  of  (2.3)  with  the  small 
scale  surface  height  spectrum  because  this  convolution  determines  [cjp^ f  ]  • 

That  is,  in  the  convolution  k  and  k  must  be  replaced  by  k  -k  and 

ox  oy  r  J  ox  x 

k  -k  where  k  and  k  are  the  new  variables  of  integration  (see  (32)  of 
oy  y  x  y  0 


If  R(0,(}))  in  (2.3)  had  been  replaced  by  its  precise  definition,  as  given 
by  (1),  the  problem  of  determining  where  kQx  and  k  appear  in  the  right 
hand  side  of  (2.3)  would  have  been  correctly  solved.  Unfortunately,  this  was 
not  done.  Instead,  R(Q,<p)  was  incorrectly  written  as  (1+CQ)  1  and,  through 
the  use  of  trigonometric  manipulations,  Cq  was  expressed  in  terms  of  k^x  and 
k^y  ,  see  equations  (24)  through  (28)  in  [1].  This  development  failed  to  rec¬ 
ognize  that  does  not,  in  general,  depend  upon  the  transform  variables  in 

(2.3).  That  is,  if  one  changes  the  transform  variables  from  k  to  k  and 

from  k  to  k  ,  C  would  still  only  depend  on  k  =  -2k  sin  8  cos  d>  and 
oy  y  o  ox  o 

kQy  =  ~2ko  sin  8  sin  <P  •  The  primary  consequence  of  this  error  is  to  provide  an 
incorrect  formula  for  R(»)  for  use  in  all  equations  following  (28)  in  [1]. 

As  will  be  shown,  this  error  fortunately  has  negligible  consequence  on  the 


numerical  results  presented  in  [1]. 


2.3  Correct  Analysis 

The  error  identified  and  explained  above  can  be  rectified  by  determining 
the  functional  dependence  of  the  right  hand  side  of  (2.1)  upon  the  large  scale 
slopes  and  since  they  are  replaced  by  the  transform  variables  kQx 

and  kQy  .  This  can  be  done  in  a  relatively  straightforward  manner  by  general¬ 
izing  Smith's  [3]  results  to  the  case  where  k^  is  at  an  angle  tt/2  —  <t>  with 
respect  to  the  y-axis  rather  than  directly  along  the  y-axis.  Such  a  general¬ 


ization  leads  to  the  following; 


p(  )^i|^£x,C£y) 


U(ctn9  -  C^cos*))  -  £Ayain<J>) 
1  +  C 


(2.4) 


where  C  is  given  by  (24)  of  [1]  and  U(*)  is  the  unit  step  function  which 
o’ 

is  one  if  the  argument  Is  positive  and  zero  if  the  argument  is  negative.  Of 
particular  note  in  (2.4)  is  the  fact  that  C^x  and  only  appear  in  the 

argument  of  the  unit  step  function.  Substituting  (2.4)  into  (2.1)  yields  the 
correct  expression  for  R(0,(j>)  ,  i.e. 


R(0,<|>) 


k  k 

U(ctn  6 - cos  <P - sin  <f>) 

D  P 

l  +  c_ 


(2.5) 


If  one  substitutes  k  -  -2k  sin  0  cos  (b  and  k  --2k  sin0sin(J>  in  (2.5) 

ox  o  oy  o 

such  as  required  in  the  determination  of  [0  °,]_  ,  then  R(0,<j>)  -  (1  +  C  )  1 

pp  U  o 

and  it  will  depend  only  on  the  angles  0  and  $  and  the  mean  square  slopes 
2  2 

and  of  the  large  scale  surface  (see  equations  (24)  and  (25)  of  [1]). 

However,  in  the  convolution  expression  for  [0^°,]^,  the  unit  step  function 
must  be  retained.  That  is,  one  must  use  the  following  relationship; 


k  -  k 
ox  X 


k  -k 

oy _ v 


k2k  cos  0  '  2k  cos  0 

o  o 


1  +  C 


(2.6) 


Equation  (2.6)  leads  to  a  completely  different  interpretation  of  the  ef¬ 
fects  of  large  scale  shadowing  on  [0^°,]^  ^rom  that  erroneously  presented  in 
(1].  According  to  (2.6),  shadowing  now  gives  rise  to  an  attenuation  factor 

(1+C  )  1  which  is  common  to  both  [o  and  [0  °f],  .  This  result  is 

o  -  PP  0  -  PP  1 

merely  a  consequence  of  the  fact  that  small  areas  on  the  surface  capable  of 
producing  strong  Bragg  scatter,  i.e.  ■  0  ,  are  as  equally  shadowed 

as  the  areas  properly  oriented  for  specular  reflection,  i.e.  “  k  /B  311,1 

=  k^/B  •  This  statement  can  be  easily  verified  by  noting  that  for  both 
sets  of  the  above  values  of  and  ,  U  (ctn©  -  cos  <{>  -  sin  (|))  =  1 


in  equation  (2.4)  and  so  R(0,<j>)  *  (1+Cq) 


Since  both  the  [o^Iq  and 


[Opp,]^  terms  suffer  the  same  attenuation  due  to  shadowing,  the  transition 
region  in  0  (where  tapp' decreases  and  [0^°, becomes  predominant) 
is  independent  of  shadowing  effects.  Since  shadowing  results  from  the  slopes 
of  the  large  scale  surface  structure  and  since  the  large  scale  slopes  are  not 
the  Important  surface  characteristic  in  determining  the  transition  region. 


this  result  demonstrates  that  the  theory  is  self  consistent. 

The  unit  step  function  in  (2.6)  serves  the  very  important  purpose  of  es¬ 


tablishing  the  limits  on  the  integrals  in  the  convolutional  expression  for 

[O  °,]  ,  see  (32)  of  [1],  However,  before  this  aspect  of  the  problem  is  con- 

pp  -L 

sidered,  it  is  worthwhile  reviewing  the  physics  behind  the  reason  for  the  unit 
step  function  in  (2.6).  The  unit  step  function  appears  in  (2.6)  because  there 
is  a  certain  range  of  surface  slope  values  for  which  the  probability  of  having 
an  incident  ray  shadowed  is  identically  one  [3].  This  result  may  be  readily 
understood  by  referring  to  Figure  2.1.  In  (a),  the  slope  of  the  surface  is 
negative  at  the  point  of  intersection  and  the  incident  ray  is  not  shadowed  in  a 
small  neighborhood  of  the  point.  In  (b),  the  surface  slope  is  equal  to  the 
slope  of  the  Incident  ray  (ctn6)  and  the  ray  is  therefore  tangent  to  the  surface 
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at  the  point  of  intersection.  In  this  case  the  point  may  or  may  not  be  shadowed 
in  a  small  neighborhood  of  the  point  depending  upon  the  surface  curvature.  In 
(c),  the  surface  slope  exceeds  the  slope  of  the  Incident  ray  and  the  ray  is 
necessarily  shadowed  by  some  portion  of  the  surface  in  a  small  neighborhood 
of  the  point.  Thus  all  points  on  the  surface  for  which  the  slope,  in  the  direc¬ 
tion  of  the  Incident  ray  projected  onto  mean  plane,  exceeds  the  slope  of  the 
incident  ray  will  be  shadowed.  Stated  another  way,  the  probability  of  such 
an  event  is  one.  This  is  the  physical  reason  for  the  unit  step  function  in 
(2.4)  and  subsequent  equations  involving  the  shadowing  function. 

In  order  to  retain  the  physical  significance  of  the  unit  step  function, 
it  is  desirable  to  deal  with  a  particular  form  of  the  equation  for  [a^°,  (0,<{»)  , 

i.e.  equation  (40)  of  [1], 


00  00 

.  „  -j  4k4(l+c  )  1  r  [ 

> . 


cos  8  £  +  k  ,2k  cos  0  £  +  k  ) 
jc  ox  o  y  oy 


U(ctn8  +£x  cos  +  £y  sin  <J>)  rpp»  (“  *»x»  ~5y) 


exp  - 


dE*  dS  -  V 


(2.7) 


where  I  is  defined  in  [1].  Equation  (2.7)  expresses  the  convolutional 
*d 

broadening  of  the  spectrum  about  the  Bragg  wavenumbers  kQx  and  kQ^  as  a 

direct  consequence  of  the  distribution  of  large  scale  slopes;  chat  is,  £x 

and  are  equivalent  to  and  .  It  should  be  pointed  out  that  the 

correct  expression  for  the  shadowing  function  has  been  used  in  (2.7).  The 

curve  in  the  £  £  -plane  separating  the  regions  where  the  step  function  is 
x  y 
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zero  and  one  is  a  straight  line,  i.e.  £x  cos  <{)  +  sin  <J>  -  -  ctn8  .  Except 

for  special  values  of  <j>  given  in  Table  I,  the  step  function  in  (2.7)  will 

consequently  give  rise  to  a  coupling  between  the  £  and  £  integrals.  More 

x  y 

specifically,  the  lower  (upper)  limit  on  the  £^  integral  in  (2.7)  is  given  by 


-  (ctn0  +  £y  sin  4>) 
cos  $ 


(2.8) 


for  cos  6  >  0  (<0)  .  From  an  analytical  point  of  view  (2.8)  represents  an 
irritating  consequence  of  shadowing.  From  a  practical  standpoint,  the  restric¬ 
tions  imposed  by  (2.8)  may  not  be  numerically  relevant  for  a  large  range  of 
incidence  angles  as  demonstrated  by  the  following  argument.  The  dominant  fac¬ 
tor  in  the  integrand  in  (2.7)  is  the  slope  dependent  Gaussian  term  which  is 
equivalent  to  the  probability  density  function  for  the  large  scale  slopes.  For 
most  practical  purposes,  the  effect  of  this  term  is  to  truncate  the  range  of 
integration  in  (2.7)  to  about  a  ± 3- sigma  excursion  from  £  -  £  =  0  .  The 


±  3-sigma  excursions  for  £  and  £  are 

x  y 


i3V/?r 


*nd 


i3*^  ■ 


respec¬ 


tively,  and  substituting  these  values  in  the  unit  step  function  argument 
yields  the  following  requirement 


±  3v  cos  <J)  ±3  sin  (J)  >  -  ctn0 


(2.9) 


for  the  unit  step  function  to  be  unity.  If  0  <  4>  <  tt/2  ,  the  "worst  case" 


situation  occurs  when  £ 


and  £ 


-3^7 


3  cos  9  +  3  vi sin  4>  <  ctn  0 


(2.10) 


The  worst  case  situation  occurs  when  the  left  hand  side  of  (2.9)  is  most  nega¬ 
tive.  Thus  for  numerical  purposes  and  all  values  of  0  satisfying  (2.10), 


one  can  replace  the  Infinite  limits  in  (2.7)  by  1  3v  £ 


and  ±  3 
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Integration 
<j>  a 


I 


ts  for  Eqn.  (2.7)  and 
tiple  of  ir/2 


ition 

upper 

Sy  Intel 
lim: 

lower 

^ration 

Its 

upper 

OO 

—00 

00 

00 

-  ctn  0 

00 

ctn  6 

—00 

00 

00 

—00 

ctn  0 
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Conditions  similar  to  (2.10)  can  be  obtained  for  other  ranges  of  <p  in  a 
straightforward  manner. 

It  is  implicitly  assumed  that  the  upper  limit  on  0  resulting  from  (2.10) 

also  satisfies  the  optical  criterion  required  of  the  large  scale  surface*, 

2  2  2 

namely  that  4kQ  cos  0  »  1  .  For  moderate  slopes  there  will  generally  be 
a  gap  between  the  maximum  value  of  0  resulting  from  (2.10)  and  the  upper 

bound  on  0  resulting  from  the  large  scale  surface  optical  criterion.  In 

«.  _  4 

this  case,  one  must  necessarily  revert  to  the  more  exact  limits  such  as  given 

by  (2.8).  What  is  happening  in  this  situation  is  that  the  ±  3-sigma  support 

of  the  slope  density  function  in  (2.7)  is  overlapping  the  region  of  the 

EE  -plane  where  the  unit  step  function  is  zero.  In  fact,  if  one  could  go  to 
x  y 

the  grazing  incidence  limit  (0-it/2),  it  is  readily  observed  from  Figure  2.2 

that  only  half  the  EE  -plane  is  encompassed  by  the  integrals  in  (2.7).  Thus, 

x  y 

in  addition  to  the  attenuation  of  [O  °,]  near  grazing  incidence  due  to  the 

r  r  ■"* 

(l+Co)1  factor,  there  is  another  reduction  factor  resulting  from  the  shadowing 
of  points  on  the  large  scale  surface  having  positive  slopes  (see  Figure  2.1). 

It  should  be  remembered  that  although  the  impact  of  the  unit  step  function  upon 
the  limits  of  the  integrals  in  (2.7)  is  somewhat  involved,  it  is  a  direct  con¬ 
sequence  of  the  rather  simple  slope-shadowing  limitations  explained  in  Figure 
2.1. 

The  analysis  presented  above  does  not  alter  the  general  composite  surface 
scattering  theory  set  forth  in  [1].  It  does,  however,  correct  and  expand  the 
theory  in  [1]  as  it  relates  to  the  effects  of  shadowing  upon  large  angle  of 
Incidence  backscattering  from  a  randomly  rough  surface.  This  additional  analy-  ' 

sis  was  necessitated  by  the  use  of  an  incorrect  functional  form  for  the  shadowing 
function  in  [1]. 

*This  criterion  has  also  been  called  the  stationary  phase  approximation  [4] 
and  the  deep  phase  modulation  condition  [5]. 
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Figure  2-2.  A  plot  of  the  regions  in  which  the  unit  step  function 
in  the  shadowing  function  is  zero  and  unity. 


2.4  Discussion  of  Numerical  Example 


The  use  of  an  incorrect  shadowing  function  in  Sections  I  through  IV  of 
[1]  was  primarily  a  sin  of  omission.  That  is,  since  nearly  all  the  results 
in  Section  I  through  IV  of  [1]  were  formal  in  nature,  it  was  not  necessary 
to  use  the  erroneous  shadowing  function  given  by  (28)  of  [1].  However,  the 
numerical  example  comprising  Section  V  of  [1]  does  require  attention  in  order 
to  properly  account  for  the  effects  of  large  scale  shadowing.  In  the  example 
presented  in  [1],  R  was  set  equal  to  unity  because  the  large  scale  slopes 
were  relatively  small.  According  to  the  analysis  presented  above,  this  step 
is  not  justified  for  all  angles  of  incidence.  The  correct  effect  of  shadowing 
is  addressed  below. 

2  2 

As  in  Section  V  of  [1],  it  is  assumed  that  T  ,  «*  T  ,(0,0)  and  the 

.  pp  pp 

2  2  2 

surface  height  spectrum  is  isotropic,  i.e.  **  / 2  and 

S(kx,ky)  -  S  (v^k^  +k^  )  .  Substituting  kcosa  ■  k^  ,  k  sin  a  -  k^  and 

dk  dk  *  kdkdct  in  (32)  of  [1]  and  using  (2.6)  above  for  the  shadowing  func- 
x  y 

tlon  yields; 


iopp-'i 


"  ht  008 


r  :,<o,o>  r  r  , 

Js -  I  J  s«0u(< 

os^ea +c)  %  y 


ctn0  +  tan  0  + 


2k  cos 
o 


os  0  / 


•  exp  j-  [4kQ2  sin2  0  +  k2  +  4kkQ  sin0  cos(ot-$)]  j 4k2  cos2  0  |  kdkda  (2.11) 

The  unit  step  function  is  unity  whenever  kcos(a-$)  >  -2kQ/sin0  .  Since  the 

minimum  value  of  cos (<x~<j>)  is  -1,  this  inequality  will  be  satisfied  for  all 

a  if  k  is  restricted  to  less  than  2kQ/sin0  .  Since  the  range  of  k  from 

2k  /sin0  to  *  for  cos(o-4>)  >  0  contributes  little  to  the  integral  in  (2.11) 
o 

2 

(because  is  small),  the  k-lntegral  limits  can  be  approximated  by 


i 
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[k^, 2k  /sin  0]  .  In  this  range  of  k  ,  the  step  function  is  unity  and  the  ot- 
integration  can  be  accomplished  with  the  following  result; 


4k2  T  2,(0,0) 

to  °t]  «  - 2_£E - ! - 

PP  1  —  2 

htcoa  0(1+Co 


2k  /sin6 
o 


—  -¥n92  )exp  {" 
+  C_)  y  \kQ  cosZ  0/  L 


(k-2kQ  sin  0) 

2~1  2  a 
4k  C04.  cos  0 
o  Sit 


k  sin  9 

.  ~~2  2 

k0^tC08  8 


(2.12) 


where  Iq(#)  i8  a  Bessel  function  of  the  second  kind.  The  right  hand  side  of 

(2.12)  may  now  be  compared  with  the  [a  °, ]  part  in  (44)  of  [1],  i.e. 

pp  x 


c.feos  0  VJ  \kt.  008^0/  L 


ht COB  6  k 


ko  h2tcos2e 


(k-2k  sin  0) 
o 


k  sin  0 


..  2  _  2  2  a  .  2  2  a  J 

4ko  CAt  cos  0  koC4t  cos  QJ 

(2.13) 


The  obvious  differences  are  the  factor  (1+Cq)  and  the  finite  upper  limit 

on  the  integral  in  (2.12).  Figure  2.3  Illustrates  how  (1+Cq)  ^  varies  with 
2 

0  for  *  0.0224  which  was  the  value  of  mean  square  slope  used  to  construct 

Figures  3  and  4  of  [1],  Of  particular  note  in  the  plot  of  (1+Cq)  ^  is  the 
fact  it  does  not  start  to  decrease  until  0  exceeds  85°;  at  87.5° 

(1+C  )-1  *  -2.5  dB.  The  value  of  0  =  87.5°  is  the  point  at  which  4k2£2  cos20  ~  10 

O  OX/ 

and,  consequently,  represents  the  approximate  limit  of  the  large  scale  theory. 
That  is,  for  0  >  87.5°  the  analysis  of  the  scattering  from  the  large  scale 
structure  on  the  surface  can  no  longer  be  accomplished  using  optical  techniques. 
The  effect  of  the  finite  upper  limit  on  the  integral  in  (2.12)  is  much  less 
dramatic.  In  fact  if  the  4kQ  cos  0  £  10  criterion  is  ignored  and  the  limit 
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ANGLE  OF  INCIDENCE  (DEG.) 

Figure  2-3.  Behavior  of  the  factor  (1+C^)"1  resulting  from 

the  corrected  shadowing  function  for  2  m  0.0224. 


of  0  “  90°  is  taken  in  (2.12),  it  can  be  shown*  that  the  finite  upper  limit 
gives  rise  to  a  3  dB  reduction  from  the  value  of  (2.13)  at  0  *  90°.  This  same 
conclusion  can  also  be  obtained  from  (2.7)  and  it  is  a  direct  consequence  of 
the  fact  that  all  positive  slopes  are  excluded  from  consideration  at  0  =  90° 

(see  Figure  2.1).  Within  the  range  of  validity  of  the  optical  criterion,  the 
finite  upper  limit  on  the  integral  in  (2.12)  gives  rise  to  much  less  attenua¬ 
tion  than  the  (1+C  )  ^  factor. 

o 

In  view  of  the  above  analysis,  it  is  concluded  that  the  numerical  results 
presented  in  Section  V  of  II]  are  correct  as  they  were  presented.  Determina¬ 
tion  of  the  spectral  division  wavenumber  is  in  no  way  altered  by  the  inclu¬ 
sion  of  the  correct  shadowing  function.  The  curves  shown  in  Figures  3  through 
6  of  [1]  are  correct  because  they  do  not  encompass  the  range  of  0  where  shadow¬ 
ing  is  important  (0~85°).  As  noted  above,  the  onset  of  shadowing  effects 

occurs  approximately  where  the  optical  criterion  (4k  cos^  0  ~  10)  is  vio- 

O  Jt 

lated,  i.e,  0  ~  87.5°  ,  for  the  numerical  example  presented  in  [1],  However, 
for  larger  slopes  shadowing  must  be  considered  since  it  will  cause  a  signifi¬ 
cant  reduction  in  app»  near  grazing  incidence. 

2 . 5  Summary 

In  the  analysis  presented  in  [1],  it  was  correctly  demonstrated  how  one 
includes  the  effects  of  large  scale  surface  feature  shadowing  on  the  backscat- 
tering  cross  sect  Lon  of  a  composite  surface  for  large  angles  of  incidence. 
Unfortunately,  an  incorrect  form  for  the  shadowing  function  was  used  in  [1] 

which  led  to  the  erroneous  evaluation  of  the  impact  of  shadowing  upon  large 
angle  of  incidence  scattering.  In  this  section,  the  correct  form  of  the  shadowing 

*To  show  this  one  can  apply  Laplace's  method  to  asymptotically  evaluate  (2.12) 
as  cos  0 -*■().  However,  it  must  be  remembered  that  the  maximum  of  the  integrand 
occurs  at  the  upper  limit  of  the  integrand  as  cos  0  -*■  0  and  this  impacts  the 
evaluation  of  the  integral  [6]. 
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function  has  been  presented  and  included  in  the  formulas  obtained  in  [1]. 
Particular  emphasis  has  been  placed  upon  the  physical  significance  of  shadow¬ 
ing  as  it  effects  large  angle  of  incidence  scattering.  It  has  been  shown  that 
shadowing  leads  to  multiplication  of  both  ]q  and  [<jppl]^  by  the  fac¬ 
tor  (1+C  )  ^  which  causes  [0  and,  thus,  0  to  go  to  zero  near 

o  PP  -t  PP 

grazing  incidence  for  a  perfectly  conducting,  randomly  rough,  composite  sur¬ 
face.  Furthermore,  there  is  another  effect  which  leads  to  an  additional  3  dfi 

h.  .  t 

attenuation  at  grazing  incidence  (6 »90‘)*.  This  effect  results  from  those 
slopes  which,  with  probability  one,  will  cause  the  point  on  the  surface  having 
these  slopes  to  be  shadowed.  At  grazing  incidence,  all  positive  slopes  are  in 
this  class. 

A  reevaluation  of  the  numerical  results  presented  in  [1]  revealed  that 
use  of  the  correct  shadowing  function  did  not  alter  any  of  the  results  rela¬ 
tive  to  the  choice  of  the  spectral  dividing  wavenumber  k^.  Furthermore,  none 
of  the  curves  presented  in  [1]  were  effected  because  they  only  encompass  the 
range  of  0  from  0  to  70°  and  the  effects  of  shadowing  were  present  for  0~85°. 
Also,  techniques  were  presented  relative  to  overcoming  some  of  the  analytical 
difficulties  resulting  from  the  use  of  the  correct  shadowing  function. 


*lt  is  reemphasized  that  exact  grazing  incidence  cannot  be  addressed  by  this 
theory  because  the  optical  criterion  assumed  of  the  large  scale  surface  fea¬ 
tures  is  violated.  The  -  3  dB  figure  is  significant  only  in  its  magnitude  rela¬ 
tive  to  the  effect  of  the  (1+C  )-^  factor. 

o 
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3.0  SHADOWING  BY  NON-GAUSS IAN  RANDOM  SURFACES 


3.1  Background 

Shadowing  of  random  surfaces  was  originally  introduced  [1]  as  an  ad  hoc 
correction  to  the  results  provided  by  physical  or  geometrical  optics  approxi¬ 
mate  theories  of  rough  surface  scattering.  Sancer  [2]  subsequently  demonstra¬ 
ted  how  shadowing  could  be  rigorously  accounted  for  in  the  optical  limit  for 
random  surfaces.  Furthermore,  he  showed  that  previously  derived  expressions 
for  the  effects  of  shadowing  based  upon  purely  geometrical  considerations  [3, 

4]  were  directly  applicable.  Using  Sancer’s  results.  Brown  [5,6]  showed  how 
shadowing  could  be  rigorously  included  in  a  formulation  for  scattering  from 
random  surfaces  characterized  by  many  scales  of  roughness,  l.e.  composite 
rough  surfaces. 

While  shadowing  theory  is  reasonably  mature,  it  has  only  been  applied  to 
jointly  Gaussian  random  surfaces.  The  Gaussian  results  are  probably  adequate 
for  the  ocean  but  they  are  questionable  for  terrain  and  completely  inadequate 
for  sea  ice  fields.  For  sea  ice,  water  first  fills  all  surface  depressions 
below  mean  sea  level  and  then  freezes.  This  elimiantes  all  surface  height 
excursions  below  mean  sea  level  and  the  probability  density  function  of  the 
surface  roughness  is  clearly  non-Gaussian.  For  these  reasons,  it  is  important 
to  extend  shadowing  theory  to  the  point  where  it  can  easily  accommodate  non- 
Gaussian  surface  statistics;  such  is  the  purpose  of  this  section. 

"l 

3.2  Analysis 

The  special  case  of  backscattering  is  chosen  to  illustrate  the  approach; 
this  minimizes  some  of  the  conceptual  details  associated  with  the  more  general 
blstatic  case.  It  turns  out  that  the  extension  of  the  results  to  the  bistatic 
geometry  can  be  accomplished  almost  by  inspection.  The  analysis  presented  by 
Smith  [4]  is  general  to  a  point  in  the  development;  however,  there  are  a  number 


of  integrations  which  must  be  accomplished  in  order  to  arrive  at  the  final 
expression  for  the  shadowing  function.  In  the  case  of  a  jointly  Gaussian 
surface  these  integrals  can  be  performed  and  a  closed-form  result  is  obtained 
for  the  shadowing  function.  For  non-Gaussian  surfaces,  the  required  integra¬ 
tions  appear  to  be,  at  best,  formidable.  The  purpose  of  this  section  is  to 

show  that  if  the  height  and  slopes  of  the  surface  are  Independent  random  vari¬ 

ables  then  the  final  expression  for  the  shadowing  function  is  drastically  sim¬ 
plified. 

For  the  reader’s  convenience  Smith’s  notation  will  be  employed  in  this 
section  and  his  Figure  1  is  essentially  repeated  here  as  our  Figure  3.1.  There 
are  three  critical  relationships  from  Smith’s  paper  [4]  which  are  required. 

If  S(F,0)  is  the  probability  that  no  part  of  the  surface  will  intersect  the 
incident  ray  (at  an  angle  0  with  respect  to  the  normal  to  the  mean  flat  sur¬ 
face)  on  its  way  to  point  F  on  the  surface  then  S(F,0)  is  given  by 


S(F,0) 


h(U  -qQ)  exp 


OO 


gCOdx 


► 


(3.1) 


where  h(»)  is  the  unit  step  function,  y  ■  ctn  0  ,  qQ  is  the  slope  of  the 
surface  in  the  y-direction  at  F  ,  and  g(x)Ax  is  the  conditional  probability 
that  the  surface  will  intersect  the  incident  ray  in  the  Interval  (T,T  +  Ax) 
given  that  it  does  not  intersect  the  ray  in  (0,t)  .  The  function  g(x)Ax  is 
determined  by  the  behavior  of  P^(£,q|F,x)  which  is  the  joint  probability  of 
the  height  £  and  y-slope  at  the  point  (x-0,y*x)  conditioned  upon  the 
height  (£q) .  and  y-slope  (qQ)  at  point  F  ;  in  particular. 


22 


gCOAx 


(3.2) 
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NORMALIZED  SLOPE,  q/w 

A  comparison  of  the  "exponential"  marginal  slope 
density  given  by  (3.13)  and  the  Gaussian. 


r 


(3.6) 


00 


Vi)P2(q)dq 


The  denominator  of  (3.5)  is  recognized  as  the  distribution  function  for  £ 
evaluated  at  £  +  VIT  ,  i.e.  F,  (£  +UT)  .  Also,  P.  (£  +UT)  in  (3.5)  is 
equal  to  the  derivative  of  the  distribution  function  evaluated  at  £  *  £q+ijt  . 
Consequently,  (3.5)  becomes  g(r)  ■  [I7F^(£q  +UT))dF^(£o  +  yx)/d(£o  +VT) 
Substituting  this  result  in  (3.1),  making  the  change  of  variable  n*£o+VT 
and  noting  that  dF^/F^  ■  d(Sn  F^)  yields 


00 

S(6,F)  -  h(ii-qQ)  exp  J  d  [5lnF1  (rj)  ]  ^ 


(3.7) 


where  Jtn  denotes  the  natural  logarithm.  Since  F^O®) 
the  following; 


1  ,  (3.7)  reduces  to 


S(0,F)  -  h(W-qo)  F1  (£q) 


(3.8) 


Remembering  that  P^(£q)  “  dF^(£Q)/d£o  ,  the  average  of  (3.8)  over  all  values 
of  £q  simplifies  to 


00 

!(PO.%>0)  -  h(y-,c)  J  {Fi<V}r/lJ  drl«o> 


(3.9) 


h(v  -  q  ) 
s(po*V6)  "  r/w  +  i 


(3.10) 


since  F^(-»)  -  0  and  T/\i  +  1^0  .  With  qQ  -  -  1/p  ,  the  shadowing  func¬ 

tion  appropriate  for  backscatter  reduces  to  the  following  simple  expression; 
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R(6) 


H 


(q/U  -l)P2(q)dq 


} 


-1 


(3.11) 


where,  in  summary,  y  -  ctn  0  ,  6  Is  the  incidence  angle,  and  ?2(q)  is 

the  probability  density  function  of  the  slopes  of  the  surface  in  the  plane  of 

incidence  defined  by  the  incident  ray  and  the  normal  to  the  mean  surface.  It 

is  interesting  to  note  from  (3.11)  that  for  normal  incidence  (0*0)  R(0)  *  1 

whereas  at  grazing  incidence  (0*Tl/2)  R(tt/2)  *  0  since  y  ■  0  and 

00 

J  q?2 (q)dq  >  0  .  Thus,  these  basic  properties  of  the  shadowing  function  are 
o 

independent  of  the  detailed  properties  of  the  slope  density  function.  One  can 
easily  verify  that  (3.11)  is  identical  to  the  results  obtained  by  Smith  for 
the  special  case  of  a  jointly  Gaussian  surface. 

The  form  of  (3.10)  compared  to  Smith's  results,  i.e.  (23)  of  [4],  suggests 
that  the  above  result  can  be  directly  translated  to  the  blstatlc  case  and,  in¬ 
deed,  this  is  the  case.  For  the  blstatlc  case,  a  generalization  of  Sancer's 
[2]  results  will  be  given.  It  should  be  noted  that  the  inequalities  involv¬ 
ing  the  angles  of  incidence  (0Q)  and  scattering  (8)  ,  just  prior  to  Sancer's 
equations  (49),  (50),  (54)  and  (55)  should  be  reversed.  With  y  * ctn  8  and 
yo*ctn8Q  ,  Sancer's  results  are  easily  generalized  by  replacing  his  Cq  by 
r(yo)/yQ  and  C2  by  T(y)/y  . 


3. 3  Example 

To  illustrate  the  above  results,  the  backscattering  shadowing  function 
R(0)  will  be  determined  for  the  exponential  joint  slope  density  function  in¬ 
troduced  by  Barrick  [7],  i.e. 

?22 (P,q)  -  exp  [-  /e(p2  +  q2)/w2]  (3.12) 
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where  w  is  the  mean  square  slope  of  the  surface  roughness.  in 

2  2  2 

(3.12)  represents  a  surface  whose  roughness  is  isotropic,  w  /2  ■  <p  >  ■  <q  >  , 
with  statistically  dependent  slopes.  That  is,  the  joint  density  function  can¬ 
not  be  expressed  as  a  product  of  the  marginal  or  individual  densities.  The 
calculation  of  the  marginal  density  P£(q)  ,  using  (3.4),  is  reasonably  straight¬ 
forward  and  the  result  is  as  follows; 

p2(q)  “  Ki lq!/w)  0.13) 

TTW 

where  K^*)  is  one  of  the  modified  Bessel  functions  of  order  one  [8].  It 
is  interesting  to  compare  this  density  with  a  Gaussian,  i.e. 

p  (q)  «  — i — exp(-q2/2w2) 

.  2  w^ 

and  this  is  done  in  Figure  3.2  where  the  normalized  densities  P£(q)w  are 
plotted  as  a  funciton  of  the  normalized  slope  q/w  .  It  should  be  noted  from 
the  plots  in  Figure  3.2  that  the  "exponential"  density  shows  a  much  greater 
probability  of  occurrence  of  small  slopes  than  the  Gaussian.  This  result  is 
in  agreement  with  one  intuitive  approach  for  generating  a  surface  character¬ 
ized  by  (3.12),  e.g.  one  strongly  filters  all  surface  height  excursions  below 
a  certain  level  to  eliminate  the  possibility  of  large  negative  height  excur¬ 
sions.  This  process  increases  the  probability  of  small  slopes  at  the  expense 
of  the  large  slopes. 

Substituting  (3.13)  in  (3.11)  and  using  tabulated  integrals  of  Bessel 
functions  given  in  [9],  the  following  closed-form  result  is  obtained  for  the 
backscattering  shadow  function  R(6)  ; 

R(9)  -  k2(x)  +  7  +  f  [Vx)Lo(x)  +  Lj_(*)K0(*)]  "j  1  0.14) 
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where  x  ■  (v^/w)ctn  0  and  the  Lq(*)  ,  n*0  or  1  ,  symbol  denotes  the  modified 
Struve  functions  [8].  Using  asymptotic  forms  for  the  special  functions  in 
(3.14),  it  may  be  readily  verified  that  R(0)  -  l(x+°°)  and  R(tt/2 )  -  0(x=0)  . 
The  modified  Struve  functions  may  be  computed  from  tables  given  In  [8]  for 
x  £  5  and  by  a  power  series  for  smaller  arguments. 

Figure  3.3  compares  (3.14)  and  the  shadowing  function  for  a  Gaussian  func¬ 
tion  obtained  by  Smith  [4]  for  a  range  of  rms  slopes.  The  shadowing  function 
for  the  exponential  joint  slope  density  is  larger  because  the  marginal  slope 
density  given  by  (3.13)  exhibits  less  likelihood  for  large  slopes  than  the  cor¬ 
responding  Gaussian  density.  That  is,  the  larger  slopes  are  the  source  of 
more  significant  shadowing. 

3.4  Summary 

The  shadowing  theory  developed  by  Smith  [4],  while  sufficiently  general 
to  deal  with  any  joint  slope  density  function,  involves  what  appears  to  be  a 
number  of  rather  complicated  integrals.  Under  the  assumption  that  the  sur¬ 
face  height  is  statistically  independent  of  the  surface  slopes,  it  is  shown 
that  Smith's  theory  can  be  reduced  to  a  single  integration  involving  the  mar¬ 
ginal  density  function  for  the  slopes  in  the  plane  of  incidence.  Using  this 
result  but  without  regard  to  the  specific  form  of  the  marginal  density  func¬ 
tion,  it  can  be  shown  that  the  backseat tering  shadowing  function  is  unity  at 
normal  incidence  and  zero  at  grazing  incidence.  Because  the  final  result  in¬ 
volves  an  Integration  or  smoothing  process,  it  is  amenable  to  the  use  of 
histogram  data  for  the  marginal  slope  density. 

This  theory  is  applied  to  an  exponential  joint  slope  density  represent¬ 
ing  an  isotropic  surface  for  which  the  slopes  are  not  statistically  indepen¬ 
dent.  The  backscattering  shadowing  function  for  the  exponential  and  Gaussian 
joing  slope  densities  are  compared  and  it  is  found  that  the  Gaussian  surface 
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ANGLE  OF  INCIDENCE,  DEGREES 

Figure  3-3.  A  comparison  of  the  backscattering  shadowing  function  for  the  exponential 
and  Gaussian  joint  slope  densities  for  a  range  of  rms  slopes  (w). 


produces  stronger  shadowing.  This  result  is  found  to  be  a  consequence  of  the 

greater  likelihood  of  large  slopes  with  the  Gaussian  density. 
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4.0  BISTATIC  SCATTERING  FROM  LOSSY  RANDOM  SURFACES 


4.1  Background 

Prior  to  the  mid-1960's,  electromagnetic  scattering  from  randomly  rough 
surfaces  was  modeled  using  either  perturbation  theory  or  physical  optics  [1], 
First  order  perturbation  theory  appeared  to  do  a  reasonable  job  of  analytically 
describing  the  scattering  process  when  the  surface  roughness  was  small  in  terms 
of  the  electromagnetic  wavelength  and  multiple  scattering  was  negligible. 
Physical  optics  produced  meaningful  results  in  and  about  the  specular  scatter¬ 
ing  direction  when  the  surface  exhibited  very  large  but  smoothly  undulating 
height  variations.  Unfortunately,  there  were  numerous  attempts  to  apply  these 
theories  to  situations  where  the  implicit  assumptions  in  the  models  were  vio¬ 
lated.  These  attempts  usually  assumed  some  surface  parameter  such  that  the 
scattering  measurements  and  the  "model”  were  brought  into  agreement.  However, 
it  was  very  quickly  recognized  that  these  attempts  were  highly  suspect  because 
of  their  failure  to  meet  certain  fundamental  principles. 

As  more  and  more  rough  surface  microwave  scattering  measurements  were 
acquired,  it  became  obvious  that  neither  first  order  perturbation  theory  nor 
physical  optics  were  individually  adequate  for  all  angles  of  incidence  and 
scattering.  Conversely,  it  appeared  that  physical  optics  seemed  to  do  a  good 
modeling  job  near  the  specular  scattering  direction  while  first  order  pertur¬ 
bation  theory  was  reasonably  accurate  for  all  other  scattering  angles.  Almost 
simultaneously,  researchers  in  the  U.S.  [2]  and  the  U.S.S.R. [3,4]  began  to 
advocate  the  combining  of  these  two  diverse  theories  in  what  was  later  to  be 
called  the  composite  surface  scattering  model.  In  this  model,  the  surface 
was  considered  to  be  made  up  of  both  large  and  small  scale  surface  features 
(height  and  spatial  wavelength)  relative  to  the  electromagnetic  wavelength, 

.  The  large  scale  surface  features  were  considered  to  be  responsible  for 
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the  physical  optics-like  scattering  near  the  specular  direction.  The  small 
scale  surface  structure  gave  rise  to  a  perturbation  field  (to  the  optical 
field)  which  was  the  dominant  scattering  mechanism  away  from  the  specular 
scattering  direction.  The  interaction  between  the  optical  and  first  order  per¬ 
turbation  fields  was  assumed  to  be  totally  dependent  upon  the  tilting  of  the 
small  surface  structure  by  the  larger  gently  undulating  features  [5]. 

More  recently,  rigorous  first  order  boundary  perturbation  theory  has  been 
applied  to  the  problem  of  backscattering  from  a  perfectly  conducting,  Gaussian 
distributed  rough  surface  [6].  The  results  of  this  analysis  indicated  that 
much  of  the  original  work  on  this  problem  could  be  rigorously  justified. 
Furthermore,  additional  insight  was  gained  in  regard  to  such  aspects  of  the 
problem  as  shadowing  (see  Section  3  of  this  report),  spectral  dichotomy,  and 
the  tilting  interpretation.  A  logical  extension  of  this  latter  theory  encom¬ 
passes  bistatic  scattering  from  a  lossy  dielectric  surface.  The  purpose  of 
this  section  is  to  present  the  details  associated  with  such  an  extension. 

Before  the  details  of  scattering  from  a  composite  dielectric  surface  are 
presented,  it  is  illuminating  to  consider  two  much  simpler  cases.  The  first 
is  backscattering  from  a  dielectric  surface  having  only  a  small  scale  rough¬ 
ness  while  the  second  case  addresses  bistatic  scattering.  The  advantages  of 
this  approach  are  that  it  leads  to  familiarity  with  the  perturbation  technique 
and  it  sets  forth  the  principles  that  will  be  used  for  the  composite  surface. 

4.2  Backscattering  From  A  Dielectric  Surface  With  Small  Scale  Roughness 

The  geometry  for  this  problem  is  shown  in  Figure  4-1.  The  mean  or  aver¬ 
age  surface  is  the  z=0  plane;  the  random  roughness  £g  superposed  upon 
this  plane  is  positive  for  £|g  > 0  and  negative  for  <  0  .  Below  the  rough 
surface  (z <  £  ),  the  relative  dielectric  constant  of  the  medium  is  e  and 
the  relative  magnetic  permeability  is  taken  to  be  the  same  as  for  free-space, 
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i.e.  y  *1  .  Above  the  rough  surface  (z > £  )  ,  the  medium  is  free-space, 
r  s 

i.e.  and  Ur  ■  1  . 

Provided  that  the  roughness  is  small  with  respect  to  the  electromagnetic 

2  2  2 

wavelength  ,  e.g.  4kQ  <<  1  where  kQ  »  2n /Xq  and  C  is  the  mean 
square  height  of  the  roughness,  the  scattered  field  Es  can  be  expressed  as 
follows; 

Es  s  6  E  +  5  E  --  (4.1) 

where  6°$  is  the  field  scattered  by  a  surface  having  no  roughness  (the  zeroth 
order  perturbation)  and  6*E  is  the  scattered  field  which  depends  on  the  rough¬ 
ness  to  first  order  only  (the  first  order  perturbation).  The  primary  assump- 

2  3 

tion  in  (4.1)  is  that  higher  order  terms  such  as  0(C  )  ,  0(£  )  ,  etc.,  are 

,  S  8 

negligible.  The  zeroth  order  perturbation  field  6°E  is  trivially  determined 
since  it  is  just  the  field  reflected  by  an  infinite,  flat  dielectric  interface. 
Both  Mitzner  [7]  and  Burrows  [8]  have  obtained  particularly  useful  expressions 
for  61E  .  The  Mitzner  result  is  more  straightforward  but  it  is  restricted  to 
small  roughness  perturbations  superposed  on  a  flat  plane.  Burrows'  solution 
for  SXE  is  somewhat  more  complicated  but  it  is  more  general  in  that  the  unper¬ 
turbed  surface  need  not  be  planar  or  even  deterministic.  For  small  scale  rough¬ 
ness  on  a  plane,  the  Burrows  formulation  requires  a  bit  more  effort  in  computing 
6XE  than  Mitzner' s  result.  For  a  composite  surface,  only  the  Burrows  result 
is  sufficiently  general  to  address  this  problem. 

At  first  glance,  the  Burrows  expression  for  o  E  appears  to  be  somewhat 
cumbersome  and  confusing.  However,  if  it  is  realized  that  the  result  is  obtained 


•Mitzner’ s  result  can  actually  be  applied  to  any  unperturbed  surface  for  which 
the  wave  equation  is  separable.  For  the  problem  considered  here,  the  surface 
is  restricted  to  a  plane. 
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from  an  application  of  reciprocity  then  the  notation  becomes  more  meaningful. 
Basically,  one  deals  with  two  incident  electric  fields  of  the  form 

”  EjL  e  \  exp(  -  j  k^ • r )  (4.2a) 

and 

Ei*'  Et'-  Eoexp(-Jk1'«r)  (4.2b) 

where  the  primed  field  may  have  a  different  polarization  and  direction  of  inci¬ 
dence  than  the  unprimed  field.  Burrows’  expression  for  the  first  order  per- 
turbation  (electric)  field  scattered  in  the  direction  -  k^'  and  polarized  in 
the  e*  direction  is  as  follows  [8]; 


61;E* 


kQ  exp (  -  jkQR) 

4uR  E  e 
o  o 


/ 


+  AB*H'  -  AH*B’  -  AD 


•E’] 


C  dS 
s  o 


(4.3) 


where  it  is  assumed  that  6  E*e'  is  measured  in  the  far-field  of  the  rough 
surface.  The  distance  R  is  measured  from  the  origin  of  the  reference  co¬ 
ordinate  system  on  the  mean  surface  to  the  point  of  observation  or  measurement 
of  fi'E^e'  and  is  the  permittivity  of  free-space.  The  fields  E'  ,  H*  , 

D  and  B  are  the  fields  on  the  unperturbed  surface  (Sq)  due  to  the  primed 

incident  field  while  AE  ,  AH  ,  AD  and  AB  are  the  discontinuities  in  the 
fields  on  the  unperturbed  surface  (Sq)  due  to  the  unprimed  incident  field. 
Thus,  to  determine  the  e^  -  polarized  component  of  the  scattered  first  order 

A 

perturbation  field  in  the  general  direction  k  one  merely  sets  e'  *e  and 

s  p 

A  _  A 

k^  ■  -ks  in  the  expression  for  the  primed  incident  field  ,  computes  the 

^  ^  ^ 

resulting  fields  E’  ,  H'  ,  D’  ,  B'  on  the  unperturbed  surface  Sq  ,  and  sub¬ 
stitutes  these  results  in  (4.3). 

For  backseat ter ing  k  ■  -k  so  according  to  the  above  receipe,  k  '  “  k. 

si  i  1 
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which  in  the  coordinates  of  the  geometry  shown  in  Figure  4-1  is  as  follows; 

k^'  •  “  kQ(  -  sin  6  cos  4>  x  -  sin  0  sin  <f>  y  -  cos  0  z)  (4.4) 

Contrary  to  previous  analyses  [1],  the  direction  of  incidence  specified  by 
.  the  angle  <|>  should  not,  at  this  point  in  the  development,  be  arbitrarily  set 

to  some  convenient  value  such  as  0  or  tt/2  .  The  reason  for  this  is  that  the 
surface  may  have  anisotropic  roughness  and  the  orientation  of  the  x  and  y 
axes  of  the  reference  coordinate  system  should  be  fixed  relative  to  this  sur¬ 
face  characteristic  and  not  the  direction  of  incidence. 

Since  the  fields  inside  the  surface  integral  in  (4.3)  are  the  fields 
induced  on  the  infinite  planar  dielectric  surface  Sq  ,  it  is  convenient  to 

further  categorize  the  problem  according  to  the  polarization  of  the  incident 
->■  •  -> 

fields.  For  both  and  E^'  horizontally  polarized,  e  and  e'  are 

/V 

orthogonal  to  the  plane  formed  by  the  unit  vectors  k^  and  n  »  z  where  n 
is  the  normal  to  the  mean  or  unperturbed  surface.  In  this  case  both  e  and 
e'  are  totally  tangential  to  the  mean  plane.  For  both  E.^  and  E^'  verti- 

A  A|  A 

cally  polarized,  e  and  e  are  parallel  to  the  plane  formed  by  k^  and 

A  _  /V 

n  m  z  • 

4.2.1  Horizontal  Polarization 

When  e  and  e'  are  tangential  to  the  mean  or  unperturbed  surface,  the 

fields  E  ,  E*  ,  D  and  D'  (when  evaluated  on  Sq)  are  entirely  tangential  to 

S  .  Since  the  tangential  component  of  the  electric  field  is  continuous  across 
i  o 

i 

an  interface,  AE-0  .  Furthermore,  since  there  is  no  change  in  across 
the  boundary  and  the  lower  medium  is  assumed  not  to  be  perfectly  conducting, 

AB  and  AH  are  both  zero  on  the  interface.  Consequently,  (4.3)  reduces  to 
the  following; 
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6  E«e *  «  - 


k  exp(  -  jk 


4ttR  E  e 

o  o 


J  (AD-1 


'K  dS 

* 3  O 


(4.5) 


The  total  E' -field  on  the  surface  Sq  due  to  E^'  is  given  by 


(1+^)  E^(Sq)  -  EQ(1+Rh)  exp  (-j£1,*L)*' 


(4.6) 


where  is  the  Fresnel  (field)  reflection  coefficient  for  horizontal  polari- 

zation  and  r^  «  xx  +  yy  or  just  r  evaluated  on  Sq.  The  discontinuity  in 
D  is  given  by 

AD  -  eQ  |  E± (z  -  0+)  -  er  ^(z- 0*)j  (4.7) 

and 

*►  »  -f 

E(z-0+)  -  (l+R^E^S  )  (4.8a) 


E(z-O-)  -  Th  E±  (Sq) 


(4.8b) 


where  T,  is  the  Fresnel  (field)  transmission  coefficient  for  horizontal 
n 

polarization.  Combining  (4.8a)  and  (4.8b)  in  (4. 7) yields 

AD  =  e  (l  +  R  -  e  T  1  E,  (S  ) 
oL  h  r  hj  i'  o' 

or 

“  -  E„[1+v£rTb]‘oex|,(-JVf.i.)8 

Multiplying  (4.7)  by  (4.9)  and  realizing  that  Th  -  1  +  yields 

AD-E*  =  -eoEo2(l  +  Rh)2(£r-l)  exp  (-j  21^-r^)  (e-e‘) 

Substituting  this  result  in  (4.5)  produces  the  desired  result  for  the  first 
order  perturbation  field  polarized  in  the  e  direction 


(4.9) 


(4.10) 
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+  k  2  E 

6lE*e'  -  -2—2- 
4irR 

where  £•§'  -  1, 


exp(  -  j  koR)(l +Rh)2(er -1)J* J  exp(  -  j  2^‘r^)  Cg  dx  dy 


(4.11) 


co 8  B-ye  -  sin2  ft 
cos0  ~  sin2  9 

and  so 

~  4(e  - 1)  cos2  0 

(1+V  (£r-D  -  . . . . -2  ^2> 

I  cos  0  +^er  -  sln20  J 

There  are  several  points  to  note  about  (4.11).  The  derivation  of  (4.11) 
was  considerably  simpler  than  the  Rayleigh-Rice  approach  [1];  this  is  because 
all  of  the  difficult  work  was  done  in  obtaining  (4.3).  Equation  (4.7)  is  an 
expression  for  the  scattered  first  order  perturbation  field,  a  more  meaning¬ 
ful  quantity  than  the  average  scattered  power  when  dealing  with  phase  sensi¬ 
tive  systems.  The  average  of  (4.11)  is  zero  because  <£  >  -  0;  however,  the 

8 

average  of  <|61E«e' |  >  or  the  Incoherent  power  is  not  zero.  At  the  beginning 
of  this  section  e'  was  specified  to  be  in  the  same  direction  as  e  ,  thus 
61E*e'  represents  6'e^  where  the  doubla-h  subscript  denotes  horizontal 
polarization  on  transmission  and  reception.  It  is  now  possible  to  examine  the 
consequences  of  cross-polarized  sampling  of  the  scattered  field.  In  this  case 
e'  is  orthogonal  to  e  ;  thus,  e  is  horizontally  polarized  and  e'  is 
vertically  polarized,  i.e. 

e  -  ■  sin  $  ft  +  cos  $ 

e'  -  e  ■  -  cos  0  cos  x  —  cos  0  sin$  y  +  sin  0  z 
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Returning  to  (4.3)  for  this  case,  it  is  noted  that  AE  is  still  zero  because 

■+  ->  ■+■ 

is  tangential  to  Sq  ,  AB  and  AH  are  still  zero  because  there  is  no 

change  in  pr  across  Sq  and  the  conductivity  is  assumed  finite,  and  the  prob- 

lem  reduces  to  evaluating  AD*E  on  SQ  .  However,  AD  has  the  direction 

e^  while  E'  is  polarized  in  the  ev  direction,  consequently,  e*e'  “e^e^O 

and  there  is  no  depolarization  by  the  surface.  This  is  just  a  confirmation  of 

the  fact  that  first  order  perturbation  theory  does  not  lead  to  a  depolarized 

scattered  field  when  the  roughness  is  small  scale. 


4.2.2  Vertical  Polarization 

The  case  of  vertical  polarization  is  a  bit  more  algebraically  involved 
because  AE  is  no  longer  zero  across  Sq  .  Both  magnetic  field  discontinu- 
ities,  AB  and  AH  ,  are  still  zero  for  the  same  reason  as  given  above. 
Thus,  (4.3)  reduces  to 


k  exp(-jkR) 
6xE*e'  -  — - 2 - 


4tt  R  E  e 
o  o 


J* [AE*D’  - 


AD  •  E’K  d  S 
s  o 


(4.13) 


For  vertical  polarization,  it  is  customary  to  use  the  incident  magnetic  field 


as  the  source.  Thus,  the  incident,  reflected  and  transmitted  magnetic 
fields  on  Sq  are  given  by 


«i  "  HoShexp(~j  V*!* 

“t  “  H0Rv  ehexp(-j 

«t  ‘  HoTvghexp(~JVr^) 


where  R  and  T  are  the  Fresnel  (field)  reflection  and  transmission  co- 
v  v 

efficients  for  vertical  polarization.  On  SQ  the  corresponding  electric  fields 
are  as  follows; 
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E 

q 


-f 

X  H 

.  q 


A 


where  q  *  i,  r  and  and  kfc  specify  the  direction  of  propagation  of  the 

incident,  reflected,  and  transmitted  fields.  Note  that  since  rx  is  on  Sq  , 


Vr-L 


k.  *r , 


which  is  merely  a  restatement  of  the  fact  that  the  angle  of  incidence  equals 
the  angle  of  reflection  and  Snell's  law  [9].  For  backscattering  and  like 
polarization  sampling  of  the  scattered  field,  the  primed  fields  are  the  same 
as  those  above. 

Although  somewhat  cumbersome  at  this  stage  of  the  development,  it  is 

desirable  to  split  the  fields  into  components  which  are  tangential  to  and 

normal  to  Sq  .  The  reason  for  introducing  this  transformation  is  that  it 

will  be  very  useful  in  the  composite  surface  development  and  it  is  therefore 

beneficial  to  obtain  some  facility  with  the  technique  on  this  easier  problem. 

The  normal  to  S  is  n  *  2  while  the  tangent  will  be  taken  as  x  =  n  x  e,  . 

o  h 

_  A  A  A 

This  particular  choice  of  x  is  convenient  because  x  •  ey  selects  the  com¬ 
ponent  of  ey  that  is  tangent  to  the  surface.  Since  AE*t  *  0  and  AD*n  *  0 

on  S  ,  (4.13)  simplifies  to 
o 


[(AE*fi)(D'*n)  -  (AD»t)(E'*t)]  ^  d Sq 

(4.14) 

where  S'  has  been  replaced  by  its  equivalent  since  like  polarization 

sampling  of  the  scattered  field  has  been  specified.  The  boundary  conditions 
AE*X  *  0  and  AD*n  ■  0  should  not  be  discarded  because  they  will  provide  some 
useful  relationships.  From  a£*t  -  0  there  results 


- 

v 


4irR  E  e 
o  o 
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r  ^  /V  A  A  1  A  V  .  A  .  A 

K  Xeh>  +  Rv(kr  X«h>]  •  T  ~  (kt  xeh)  .  T 


(4.15) 


while  A5»fi  “  0  yields 

[(ktxah)  +  Rv(krxSh>]  -  /r  Tv(ktx^h).  S  (4.16) 

where  AE  -  E.  +E  -  E  and  AD-D.  +D  - D  have  been  used  along  with  (4.15) 

1  r  t  l  r  t  »- . 

and  D  *  6  E  • 

It  is  now  necessary  to  determine  the  field  quantities  inside  the  integra- 
tion  in  *4.14).  The  quantity  AE*n  can  be  reduced  to  the  following  form 
through  the  use  of  (4.15); 


AE*n  -  — ml -~r~  H 
1  eo  ° 


(er  -1)  exp  (  -  j  k1*rx)(kt  xeh)*n  (4.17) 


while  (4.16)  simplifies  D' *n  to 


D -n  -  -  /por  Hq  Ty  /T  exp  (  -  j  k±4x  )  (&t  x \)  *6 


(4.18) 


so 


(AE 


■n)  (D'  *n)  -  UQ  Hq2  Ty  (er -1)  exp  (  -  j  2^*?^)  [(^t  x  ^h>  *n]2  (4.19) 


Through  similar  manipulations. 


-f  A  A  A 


a5*t  -  -  /yQ£o  Hq  — —  (1  -e^)  exp  (  -  j  k^r^)  (kt  xeh)*T 

/z~ 

r 


(4.20) 


and 


*  a  K  \  ■+  +  ,*• 

Hq  —  exp(- jk^r^Ckj. 


o 
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x  ah).T 


(4.21) 


so 


T 

(AS»t)(E'*t)  -  uoHq2  ~~  (1  -er)  exp  (  -  j  lc^r^)  [(kt  *eh)*T]  (4.22) 


Combining  (4.19)  and  (4.22)  and  completing  the  unit  vector  operations  yields 


(AE*n)  (D'  *n)  -  (AD»t)  (E1  •!) 


2  2  (er-^  r  2  i  -*■ 

v.  v  -^-[et+(v1)8ln  e]e*pl-321irrj.) 

r  ...  (4.23) 


where 


2e  cos  6 

T  “  - - - 

v  j - ~ 

e  cos  6  +  /e  -  sin  0 
r  r 

Substituting  this  result  in  (4.14)  and  recognizing  that 

and  e  --k.  ,  the  final  result  is  obtained 

v  i  h 


E 

o 


H 

o 


HoTv  — ^~2~  [er  +  (er-1)8in2^jJexp(-j2k1»rL)Csdxdy 

(4.24) 

Essentially  the  same  remarks  apply  to  the  vertically  polarized  scattered 
field  as  for  the  horizontal  case.  In  addition,  it  should  be  noted  that  if 
(4.11)  and  (4.24)  are  converted  to  0°  or  the  scattering  cross  section  per 
unit  area  according  to 


« u  ■  k°e*p;:^) 


£ 


a°(0,4»)  *  lim  lim 
R-*»  A -*• 


f  4tt  R2  <  1 6 1 E 1 2> 

1  A  E2  i 


where  A  is  the  illuminated  area,  the  result  is  identical  to  the  result 
obtained  by  Peake  using  the  Rayleigh-Rice  approach  [1], 
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4.3  Bistatic  Scattering  From  A  Dielectric  Surface  With  Small  Roughness 

For  bistatic  scattering  the  unit  vectors  specifying  the  directions  of 
incidence  of  the  unprimed  and  primed  fields  are  given  by  (see  Figure  4-1) 

A  A  A  A 

k^  ■  -  sin  9^  cos  <J>^  x  -  sin  9^  sin  y  -  cos  9^  z 

(4.25) 

kj  ■  -  sin  9  cos  x  -  sin  9  sin  4>  y  -  cos  9  z 

i  s  s  8  s  s 


•>  A  ^  _ 

and  k.  ■=  k  k.  ,  k  *  k  k,  .  The  unit  vectors  specifying  the  directions  of 
i  oil  o  i  r 

horizontal  and  vertical  polarizations  for  the  unprimed  and  unprimed  fields 
are  as  follows; 


e^  “  -  sin9^  fc  +  cosij)^  $ 

e '  **  -  siiwj)  x  +  cos<|»  y 

h  S  8 

(4.26) 

ey  *  -  cosS.^  cos4»i  x  -  cosS^  sin^  y  +  sinO^  z 

e'  -  -cos9  cos*  x  -  cos9  sin*  y  +  sin9  z 

v  s  s  a  s  s 

The  normal  to  the  unperturbed  surface  is  n  *  z  while  the  tangents  to  the 
surface  for  vertical  polarizations  are  given  by 


A 

T 


n  x  e. 


A  A 

n  x  eh' 


(4.27) 


Equations  (4.25)-(4.27)  describe  the  basic  quantities  that  will  be  required  in 
this  section. 


4.3.1  Horizontal  Polarization 

Since  AB  and  AH  are  zero,  because  there  is  no  change  in  ur  across 

S  and  the  conductivity  is  finite,  (4.3)  becomes 
o 


exp  (-  j  kQ  R) 

4tt  RE  e 
o  O 


I 


[AE*D'  -  AD»E* ]  C  d  S 
s  o 


(4.28) 
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Furthermore,  AE  is  also  zero  across  Sq  because  E  is  tangential  to  Sq  . 
Thus,  (4.28)  becomes 


5'K’ 


k  exp (  -  j  k 


4lTR  E  e 
o  o 


-! 


[AD*E']  C  dS 
8  0 


(4.29) 


On  the  surface  S  ,  the  incident,  reflected 
o 

(unprimed  and  primed)  are  as  follows; 

\  m  I0.W<-lk1-r_)«h 

Er-EoRh«p(-^r-tL)eh 


\  -  EoThe,p(-jkt.rj_)6h 


where  also  on  the  surface  Sq  , 


k  *r.  ■  k  *r  -  k  *r 
i  ±  r  j.  t  x 


,  and  transmitted  electric  fields 

“  Eoexp(-j^.^)Sh' 

I;  -  EoRh'exp(-j^.^)gh' 

Et'  -  Eq  Th'  exp  (  -  j  kj  •^)eh’ 

(4.30) 

k,’  -  k  '  -r  -  k.'  -r  (4.31) 

i  O.  r  X  t  x 


and  the  same  notation  as  introduced  earlier  has  been  continued.  The  fields 

AD  and  E'  on  S  are  as  follows; 
o 

^  •  eo  V1  +  V  er  Th)  exp  (  '  J  W 

E'  -  Eq(1  +Rh' )  exp  (-  j  k^r^)  eh' 

Using  the  fact  that  1  +  ■  T^  ,  the  product  AD»E'  becomes 

AD*E'  -  eQ  E*  (1  +R|i)(l  +Rh')(l-  er)(gh’gh  )  exp  [“j  ^i  +^i  * 


and  substituting  this  result  in  (4.29)  yields 


k  exp (-jk 

-  -S  ° 


4ttR 


—  VMVa"b>(e  rl)<Sh,8h)J  J  e«p[-j(Jl4p.rJ;ild.4y 


(4.32) 
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This  result  can  be  easily  translated  into  the  angles  (8.  »<{>,)  and  (6  ,<p  ) 

11  «  S  8 

by  the  use  of  (4.25)  and 


W  "  cos(4>±-  4>s) 


along  with 


2  cos  01 


2cos  0 


1  + 


cos  0, 


+  A  -  sin2  0. 


1  +  Rj 

n 


s 


CO  8  0 


sin2  0 
r  s 


For  ease  of  comparison,  it  should  be  noted  that  most  results  similar  to  (4.32) 

A 

express  the  direction  of  scattering  as  k  which,  in  the  above  notation,  is 

s 

A  J 

-  .  Equation  (4.32)  yields  6  E  .  For  cross  polarized  sampling  of  the 

scattered  field,  the  problem  becomes  somewhat  more  involved  and  it  will  be 
discussed  in  Section  4.3-3.  It  should  be  noted  that  when  <{>g  =  (backscat- 
tering) ,  (4.32)  reduces  to  the  result  obtained  in  Section  4.2.1. 


4.3.2  Vertical  Polarization 

Since  n*E  and  t*D  are  discontinuous  across  the  unperturbed  surface, 
(4.3)  reduces  to 


ko2exp(-jkQR) 

4irRE  e 
o  o 


(AD*t)  (E*  •t’  )  (t*t'  )J  £s  d  Sq 


(4.33) 


A  A.  A 

where  the  scalar  product  T*T  must  be  included  because  the  unit  vectors  T 
and  T*  are  not  necessarily  parallel,  e.g.  see  (4.27).  The  incident,  reflec¬ 
ted,  and  transmitted  unprimed  and  primed  magnetic  field  quantities  are  given  by 


n 

■H 

Hoe,p(-jt1-?x)Sh 

5i 

-  H0€xp(-Jki 

-► 

H  - 
r 

H0Rve*p(-j£t-\)4h 

S' 

r 

-  H0Rv'exp(-3k;.\)eh' 

,  *>■  k  . 

■+■  . 

Ht- 

Ht 

*  H0Tv’exp(-jkt'*r_L)eh' 

(4.34) 
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while  the  corresponding  electric  fields  are 


AE 


"  £  o  r—  r 


JE~ 

'  V 


l)exp(  -  j  k4*r  )  (kfc  x  eh>*n 


(4.38) 


and 


Ho  /T  Tv'  exp  (  -  j  k^-rj  (kj  x  eh*  )  .fi  (4.39) 


In  a  similar  fashion,  (4.36)  is  used  to  simplify  the  expressions  for  AD»x 


A 

and  E  »x  with  the  following  result; 


=  ^M0  Hq  eXp(-j  k14x)(ktxgh).^ 

/e 

r 


(4.40) 


and 


E' 


A  fir  t  • 

•T’  Hq  exp(-j  xgh').^  (4.41) 

r 


Substituting  (4-38)  -/(4.41)  in  (4.33),  noting  that  x»x'  ■  cos(<J>  ~<t>  )  , 

1  s 

and  simplifying  the  unit  vector  operations  yields  the  following  result  for 


'  ; 


5‘E‘e  ' 


ft-  H 
»  e  o 


exp(-jk  R)  (£  -1) 

- —  T  T  »  — 


e  o  4irR  v  v  2 
o  e 

r 


<  £  sinS^  si 

ir  1 


in6  +  /{z  -sin^Q . )  (e  -sin^B  ) 


2  2 

3.  )(e  -  • 

r  l  v  r  s 


•  COs(<J>^  -  <j>a  )j  J  J  exp  [-  j  (k±  +ki' )  *rx)3cs  dxdy 


(4.42) 


where  E  »  H  /y  /e  and 
o  o  o  o 


T  =  1  +  R  = 

V  V 


2  c  cosB . 
r  i 


e  cos  0 , 
r  i  r 


+  l/c  -  sin^0. 
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and  T'  -  T  (6, -*-0)  .  Equation  (4.42)  is  6>E  ,  and  it  is  easily  shown 

to  reduce  to  the  result  for  backscattering  in  Section  4.2.2.  when  d>  *  d> 

i  s 

and  0^  -  0g  .  Once  again  it  should  be  emphasized  that  the  most  difficult 

part  of  obtaining  (4.42)  is  evaluating  the  terms  (k  x e,  )*x  and  (k  '  x  g  ')• 

t  h  t  h 

in  (4.40)  and  (4.41).  As  noted  previously  this  is  a  consequence  of  the  fact 
that  the  difficult  analysis  was  finished  once  (4.3)  was  derived  and  the  actual 
evaluation  of  (4.3)  is  very  straightforward.  Finally,  comparing  the  a° 
values  resulting  from  (4.32)  and  (4.42)  with  the  corresponding  results  obtained 
from  the  Rayleigh-Rice  theory  [10]  shows  complete  agreement. 

4.3.3  Cross  Polarization 

As  shown  in  Sections  4.2.1  and  4.2.2,  depolarization  for  scattering  by 

small  scale  roughness  is  a  second  order  effect  in  the  plane  of  incidence. 

If,  however,  the  scattered  field  outside  of  the  plane  of  incidence  is  computed, 

it  will  be  found  to  have  a  nonzero  cross  polarized  component.  This  result  is 

simply  a  consequence  of  the  fact  that  the  unit  vectors  e^'  and  ej  are 

not  fixed  with  respect  to  the  surface-centered  coordinate  system  and  they 

change  their  directions  as  the  observation  point  moves  out  of  the  plane  of 

incidence.  This,  of  course,  is  a  purely  geometrical  effect  and  it  has  nothing 

to  do  with  any  change  in  the  basic  scattering  mechanism. 

The  derivation  of  the  results  follows  essentially  the  same  pattern  as 

set  forth  in  the  previous  sections.  There  is  one  point  that  should  be  noted 

because  it  simplifies  the  algebra  somewhat.  For  the  case  of  the  incident 

field  horizontally  polarized  (@^)  and  the  scattered  field  vertically  polarized 

(g^' )  ,  the  unprimed  field  quantities  should  be  obtained  from 

E .  “  E  exp(-j  k.  »r,  )e.  while  the  primed  fields  should  be  obtained  from 
i  o  l  X  n 

H^'  ■  Hq  exp (  -  j  k^*r^)e^  .  For  the  incident  field  vertically  polarized  (ey) 
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<  w 


and  the  scattered  field  horizontally  polarized  (e^' )  ,  the  unprimed  fields  are 
obtained  from  «  Hq  exp (  -  j  k^*r^)e^  while  the  primed  fields  are  to  be  de¬ 

rived  from  E^'  =  EQexp(-j  kj*r^)§h'  •  This  approach  is  consistent  with  the 

technique  of  obtaining  all  field  quantities  from  the  horizontal  (e,  or  e  ' ) 

n  h 

field  for  planar  surface  reflection. 

With  E  =  H  /y  /e  ,  the  following  expressions  for  61E,  .  and  61  E  ,  , 

O  O  O  O  or  Tiv  vh 

result; 


<5^ 


hv' 


61E»e' 


kQexp(-jkoR) 
ir  R 


V 


cos0.  cos0  sin(<j)  .  )  (e  -1  )v4  -sin^0 
i  s  s  i  v  r  r  s 


(cos0  +  \/z  +sin^0.  )  (  e  cos  0  +  Vc  +sinE0  ) 
i  r  l  r  s  r  s 


If 


exp  [-  j  (k1+k1’  )  *rx]  t;8  dx  dy 


(4.43) 


6  Evh* 


k  exp (- 


^’4 

ttR  le  o  I 


COS0.COS0  sin(c|>  -<b . )  (e  -1)  l/e  -sin^0, 
i  s  ri'  r  r  i 


(cos0  +^4  -sin^0  ) (e  cos0  -sin^0 . ) 
sr  sr  ir  i. 


•  J J  ex  p  [-  j  (k£  +  k^  )  •  rj  Cs  dx  dy 


(4.44) 


where,  in  summary,  the  angles  are  defined  in  Figure  4.2,  r^  =  xx  +  yy  ,  and 

*>■  “►  *  A  A 

and  k^  are  defined  as  k^k^  and  kQk^  ,  respectively,  where  k^  and 

/v 

k^'  are  given  in  (4.25). 

When  comparing  (4.43)  and  (4.44)  with  the  cross  polarized  scattered  fields 

resulting  from  the  Rayleigh-Rice  approach  [10],  should  be  set  equal  to  n  . 

The  expression  for  6 1 ,  agrees  with  the  results  in  [10,  pg.  706].  The 

expression  for  6JE,  ,  is,  however,  the  negative  of  the  a,  coefficient  in 
nv  nv 

[10,  pg.  706,  eqn.  9.1-69].  Normally  this  difference  is  not  important  because 
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61 ia  squared  and  then  averaged  to  find  the  incoherent  power.  However, 

if  one  is  dealing  with  circular  polarization  the  sign  does  become  critical. 

To  resolve  this  issue,  a  special  case  can  be  constructed  whereby  ot^y  should 

agree  with  or  S1  ■  61  E^  .  This  special  case  involves  taking 

d>,  ■  it  and  4=0  and  8  **  0  in  the  expression  for  a  or  61  E  ,  and 

Ti  rs  s  w  w 

comparing  this  result  with  a^v  and  61  E^v,  for  <j)^  =  u  ,  <j>g  =  tt/2  ,  and 

0^  *  0  .  In  this  special  case  both  61Eyy,  and  61E^v,  should  be  polarized 

in  the  -  x-direction  .  Comparing  (4.42)  and  (4.43)  for  this  special  case  shows 

that  indeed  5'e  ,  *  <5*  E.  ,  .  However,  evaluating  a  and  a  from  [10] 

w  hv  w  hv 

results  in  *  -a^y  ;  consequently,  there  does  appear  to  be  a  sign  error 

in  the  expression  for  a  and  (4.43)  is  correct. 

This  section  completes  the  development  for  scattering  from  a  dielectric 
surface  having  only  a  small  scale  roughness.  Once  again  it  should  be  emphasized 
that  the  purposes  of  Sections  4.2  and  4.3  are  (1)  to  check  the  Burrows  pertur¬ 
bation  approach  against  the  conventional  Rayleigh-Rice  results  and  (2)  to  il¬ 
lustrate  the  actual  mechanics  of  evaluating  the  Burrows  expression  for  the 
first  order  perturbation  field.  Hopefully,  this  latter  purpose,  if  achieved, 
should  considerably  simplify  the  transition  to  the  composite  surface  case. 

4,4  Bistatic  Scattering  From  A  Dielectric  Surface  With  Composite  Roughness 
For  small  scale  roughness  superposed  on  a  planar  surface,  the  Burrows 
perturbation  formula  (4.3)  is  particularly  easy  to  evaluate.  This  results  from 
the  fact  that  one  deals  with  the  fields  on  an  infinite  planar  surface  and,  for 
such  a  surface,  the  fields  are  easily  described  and  related  through  the  Fresnel 
coefficients,  Snell's  law,  and  the  equality  of  the  angles  of  incidence  and 
reflection.  For  a  composite  surface,  the  unperturbed  surface  is  not  planar 
but  it  is  assumed  to  be  very  gently  undulating.  More  specifically,  the  unper¬ 
turbed  surface  is  actually  defined  such  that  it  contains  no  spatial  frequency 
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components  which  are  smaller  than  Y^0  ,  where  Y  is  a  constant  which  is 
greater  than  unity.  Of  course,  it  is  desirable  to  have  y  as  large  as  pos¬ 
sible  but  this  is  not  always  practical  since  the  small  scale  height  must  satisfy 
2  2 

4k  C  «  1  [6]  .  However,  if  y  can  be  made  sufficiently  large  then  the 

o  s 

scattering  from  the  unperturbed  surface  can  be  treated  using  physical  optics. 

Physical  optics  assumes  that  the  surface  may  be  considered  to  be  locally  planar 

and  the  fields  on  the  surface  can  be  accurately  approximated  using  Fresnel 

theory.  This  approach  is  recognized  to  be  essentially  the  same  as  the  small 

scale  roughness  on  a  planar  surface  problem.  The  one  important  difference  is 

that  for  the  gently  undulating  unperturbed  surface,  the  local  normal  is  no 

longer  entirely  z-directed  and,  in  fact,  depends  upon  the  slopes  of  the  large 

scale  surface.  This  means  that  one  must  construct  a  local  coordinate  system 

0 

on  the  undulating  unperturbed  surface  and  compute  the  surface  fields  required 
in  (4.3)  in  terms  of  this  system.  This  must  be  done  for  both  the  unprimed  and 
primed  fields  because  they  have  different  angles  and  directions  of  incidence 
for  the  general  bistatic  case. 

For  the  unprimed  fields,  the  important  unit  vectors  are  and  n^ 

which  is  the  normal  to  the  large  scale  or  unperturbed  surface.  These  two  quan¬ 
tities  are  important  because  they  form  the  local  plane  of  incidence.  One  next 
constructs  unit  vectors  and  e^  which  are  orthogonal  and  parallel, 

respectively,  to  and  n^  .  These  unit  vectors  are  also  horizontally  and 

vertically  polarized,  respectively,  with  respect  to  the  local  plane  of  inci¬ 
dence.  Any  arbitrarily  polarized  unprimed  incident  field  can  now  be  decomposed 
into  components  parallel  to  e^h  and  e^  since  the  incident  field  must  be 
transverse  to  .  The  unprimed  field  quantities  required  in  (4.3)  can  then 
be  computed  as  in  the  previous  sections.  The  exact  same  construction  of  e^ 
and  and  the  decomposition  of  the  primed  incident  field  must  be  performed 
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in  order  to  compute  the  primed  fields  required  in  (4.3).  Fortunately,  this 
is  easily  accomplished  by  simply  changing  6.  to  8  and  <j>.  to  $  in 
the  unprimed  quantities. 

Before  getting  into  the  actual  details,  there  are  a  few  other  points  that 
should  be  noted.  All  of  the  above  noted  manipulations  are  goint  to  lead  to 
the  following  changes  in  6*E  obtained  in  equation  (16)  of  [6],  First,  the 
factor  Tppl  is  going  to  depend  on  ,  the  angles  and  directions  of  inci¬ 
dence  and  scattering  (8^  ,  ,  0^  ,  ,  and  the  slopes  of  the  large  scale  or 

unperturbed  surface  ’  T^e  on^y  other  change  is  that  the  exponen¬ 

tial  inside  the  surface  integral  will  become  exp [  -  j (k^ + k' )  • r^]  where 
r^  =  xx  +  yy  +  ?£  z  because  of  the  generalization  to  bistatic  scattering. 
Except  for  correcting  [6]  to  properly  include  shadowing,  as  detailed  in  Section 
2,  all  other  aspects  of  the  solution  presented  in  [6]  remain  the  same.  Com¬ 
bining  61E  from  this  analysis  with  Sancer's  result  [11]  for  essentially  S°E 
yields  the  total  scattered  field.  Furthermore,  it  should  be  expected  that  for 
backscattering  the  dielectric  nature  of  the  surface  should  have  an  almost 
negligible  effect  upon  the  wavenumber  at  which  the  surface  height  spectrum 
is  partitioned  into  large  and  small  scale  sub-spectra.  Finally,  because  there 
is  no  discontinuity  in  magnetic  properties  across  the  unperturbed  surface  and 
the  conductivity  is  assumed  to  be  finite,  Ab  and  Ah  in  (4.3)  will  be  zero. 
This  fact  holds  true  regardless  of  any  tilting  of  the  locally  planar  surface. 

The  first  task  at  hand  is  to  construct  e^  and  e^v  .  Since  e^  is 

A  A 

orthogonal  to  both  k^  and  n^  ,  it  is  given  by 


ki*fi£ 

l*i*\l 


(4.45) 


where  (t.  •  k.  x  +  k,  9  +  k ,  z  and  the  k.  ,  q  ■  x,y  and  z  ,  are  obvious  from 
i  ix  iy  iz  iq  H  J  * 

(4.25).  The  normal  to  the  large  scale  surface  is  given  by 
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(4.46) 


-Six*  ~  +  * 

Vi  +  c0„2  +  5.  2 


jj  ■  A  -  A 

"nAxX  +  n£yy  +  n£zZ 


’£x  •  *£y 

Expanding  Che  cross  product  in  (4.45)  yields 


S*h  ‘  hHx  8  +  hiy  ^  *  hl2  2 


where 


(4.47) 


(k.  n0  -  k.  n.  ) 
h  _  iy  Iz  iz  ly 

l£i  X  *1* 


lx 


% 


(kiz  n£x 


kix  nlz) 


l£i  **l\ 


(k.  n„  -  k.  n0  ) 

h _ i*-M _ i 2_li_ 

u  l*i  *  «»l 

and  e^h  is  completely'  determined.  For  the  unit  vector  ,  things  are 

a  bit  more  involved  because  it  must  be  in  the  plane  formed  by  k  and  n. 

X  Xi 

and  also  orthogonal  to  .  For  the  reflected  and  transmitted  fields  k^ 

A  A 

goes  to  k  and  k  .  This  will  not  change  the  direction  of  e0,  because 
r  t  xn 

k^  ,  and  are  all  coplanar.  This  will,  however,  alter  the  direction  of 
e^v  ;  this  is  easily  understood  by  noting  that  k^  ,  e^  ,  and  e^v  form  a 
mutually  orthogonal  triad  of  unit  vectors.  Thus,  if  does  not  change 

direction  but  k^  does  then  e^v  must  necessarily  change  direction.  What 

.  /N 

this  means  is  that  we  must  find  a  new  en  for  each  value  of  k  .  This  is 

xv  q 

easily  done  by  the  following  equality; 


“  @£h  *  kq 


(4.48) 


where  q  -  i,  r,  and  t  .  Note  that  it  is  not  necessary  to  divide  the  rhs  of 

A  A 

(4.48)  by  the  magnitude  because  it  is  unity,  i.e.  e^  and  k^  are  mutually 
orthogonal  by  (4.45).  The  unit  vector  8^  may  also  be  written  as 
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(4.49) 


*  q  q  n  *  n  />» 

Bo  ■  v.  x  +  v„  y  +  v,  z 
xv  Ax  Ay  Iz 


where 


•(k  h.  -  k  h0  ) 
qy  J Iz  qz  ly 


-(k  h0  -  k  h„  ) 
qx  £z  qz  lx 


v0q  -  -(k  h.  -  k  h#  ) 

Iz  qx  ly  qy  lx 

and  e^q  is  completely  determined  for  q  *  i  ,  r  ,  and  t  .  For  the  primed 
fields,  ef',  and  e«  are  obtained  by  merely  replacing  k  by  k  '  ,  q  *=  i, 


q  q 


r,  and  t,  in  the  expressions  for  e^h  and  e ^  . 


The  incident,  reflected,  and  transmitted  unprimed  fields  on  the  unper¬ 
turbed  surface  will  now  be  decomposed  into  locally  horizontal  and  vertical 
components.  If  the  incident  field  is  of  the  form  E^  =  exp  (  -  j  k^ * tj)  e^ 

on  the  unperturbed  surface  where  e  is  its  polarization  direction  then  E  , 

d  1 

-f 

E  and  E^  can  be  written  as  follows: 
r  t 

I  »  e5  i0.  +  E1  8* 
i  Ti  Ah  v  Iv 


E  -  E5  8.  +  Er  8.r 

r  h  £h  v  £v 


(4.50) 


where 


•*  _t  ^  .  _t  *  t 
Et  ‘  *h  eHh  +  EV  H, 


(4.51) 


and  q  »  i,  r,  t  .  The  corresponding  magnetic  fields  are  given  by 


H  -  t/S"  k 
q  Til.  i 


e  ^ 

—  k  x  E  — - 
^o  q  q  «qt 


(4.52) 


56 


where  5  -  ^  for  q  -  i  and  r,  and  6tfc  -  1  .  Expanding  (4.52)  yields 


vV§  (- stv  ♦  ^  ) 


From  Fresnel  theory,  the  e^-component  of  Hr  is  equal  to  R  "times  the 

/v  +  V  \ 

e^-component  of  ,  so 


Er  -  R  E1 
v  v 


(4.53) 


Similarly,  the  e^-component  of  Ht  is  equal  to  Tv  times  the  e^-component 

-  \ 

of  ,  so 


E*  ■  — —  E^ 

v  /Tx  v 


(4.54) 


Equations  (4.53)  and  (4.54)  can  now  be  used  in  (4.50)  to  express  all  of  the 
fields  in  terms  of  E*  and  E^  ,  i.e. 

h  ■  <  •»  +  <  £tv 


►  i  A  i  a  r 

l  -  R,  ET  e..  +  R  E  e0 
r  n,.  h  Ah  vj,  v 


f* '  *» + #  <  V 

r 


(4.55) 


•The  subscripts  "Jt"  on  Ry,  Rh,  Tv,  Th  means  that  the  angles  in  the  appropriate 
Fresnel  formulas  must  be  defined  with  respect  to  the  normal  to  the  large  scale 
surface,  fia  . 

*>  *7 


Id. 

where  E^  end  Ey  are  given  by  (4.52)  with  q»i  .  The  corresponding  D-fields 

•V  4-  4 

are;  D.  “  e  E,  ,  D  ■  e  E  ,  D  ■  e  e  E  .  The  primed  fields  may  be  obtained 
i  oi  r  or  t,  o  r  t 

from  (4.55)  and  (4.51)  by  replacing  k^  by  k^'  ,  e^  by  e^»  by  e^J*  and 
by  changing  e  to  whatever  scattered  field  polarization  is  to  be  sampled, 
say,  efe  . 

The  appropriate  form  of  (4.3)  is 


kQ2  exp(  -  j  kQR) 

4ttR  E  e 
o  o 


AE-fiA) (D* 


S 

o 


-  AD 


*] 


C  d  S 
s  o 


(4.56) 


where  the  shadowing  factor  has  been  temporarily  omitted  from  the  integrand 
since  it  can  be  added  at  the  end  of  the  development.  The  (AE»n^)  is  equal 

*4  *4*  ^ 

to  (E^ +Er  -  Et)*n^  or  using  (4.55)  and  noting  that  n^e^  *  0  » 


+  R 


*  r  * 

Vj^  eAv  nl 


A  t 
e£v 


•n 


*) 


(4.57) 


(4.57)  can  be  simplified  somewhat  by  using  the  relationship  that  results  from 
AD*fi  =*  0  .  The  final  result  is 


(g  -1) 

AE*n  -  E1  — ^ -  T  (e^-n.) 

i.  v  r~—  v  lv  r 

«- 


(4.58) 


For  15'  »n^  , 


,1' 


D'*n0  -  E*  e„/E7  t;  (e^'*nB) 

Z 


In  —  u  W  X  VO  A  -Ufl  J 

Z  v  o  r  Zv  Z 


(4.59) 


so  the  product  of  (4.58)  and  (4.59)  can  be  written  as  follows 

-  E02.xp[-J(£14;).ri]e0fcr-1)T'  T  V  <*£  •*,)  «-60) 

X  X* 


>  -- ** •  *• 
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It  can  be  shown  that 


lv  **■  /F 


elv  *n£  " 


ain6U 


where 


sin0Jll.  Vl-Ck.-fi/ 


sin  6^  - 


1 1  -  (k^»n^)2 


so  (4.60)  becomes 


-  (C  •  ^  -  (£  '  •  n  1  ^ 

-  (^7-  .4- 


(4.61) 


For  the  remaining  term  in  (4.56),  the  Important  parts  are  (AD)p  811(1 

(E')  where  the  p-subscript  denotes  tangential  to  the  large  scale  surface, 

■>  A 

l.e.  the  normal  component  of  D  is  continuous  across  the  boundary  so  AD*n^*0  . 
The  tangential  components  of  the  AD-field  can  be  found  by  decomposing  Ad  into 
components  directed  along  e^  and  ,  where  »  n^  *  e^  ,  i.e. 


(^)p  -  <4S.alh)  8th  +  (40-^)  T* 


(4.62) 


Using 


48  *  £o  +  eoEv1(SH1v+S£tv"Tv’^  S*) 
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and  simplifying  this  expression  with  the  aid  of 


i  V 


+  R 


,A  r  a 

vt  <•«»•  V 


&  /A  t 

^  <e^T) 


which  results  from  (AE)  ■  0  yields 

P 


,‘t),-£.‘»i<uv(er-listt 


— 2-  E  1  T  (e  - 1)  (e-1  *T) 

/T"  v  vf,  ,  tv 

(4.63) 


The  appropriate  expression  for  (E')  is 


<l 


■>p  ■  *b1’  t„;  zu> + \ 


!• 


T  ' 

V«,  -  t* 

7F  ^ 

r 


Ti)Ti 


(4.64) 


since  on  the  unperturbed  surface  E'  -  »  Ej.  .  The  relationship  for 

(AD)  (E*)  is  obtained  by  taking  the  dot  product  of  (4.63)  with  (4.64).  Com- 
P  P 

bin in g  this  result  with  (4.61)  and  substituting  into  (4.56)  yields  the  follow¬ 
ing  result; 


6*E 


ab 


Eoko  exp(-jkR) 
irR 


f  fT ab(hX’hy)I(x>y)exp  ["■ i  (Vai  >  Cs  &  dy 


(4.65) 


where  6*E  ^  is  the  scattered  first  order  perturbation  field  for  an  incident 
polarization  e^  and  scattered  polarization  e^  ,  I(x,y)  is  unity  on  the 
illuminated  parts  of  the  surface  and  zero  for  the  shadowed  parts,  and 


JItSi 
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(£  -1)  |  2  T 

4  V+h 

x  y 


(Sa‘^v>  (gb  -aiV  >  (  -pr)  er  sln  0U  8ln  6U 


+  /£r-8ineU  /er’8in0ii 


27”  ^7-2.,  f  - 


ain  0£i  sin  0^ 


+  (vev(vevv>i 


w +  r  *±w*  — rr 

sin  0£±  sin  0j^  (ea’eJlv)  (eb*eAh)  Th„  (  e_  J  ^r  "  8in  0W 


ll\  r 


./V  A  1#  -  -A  A  j  ^ 

+  (eb*e£v)(Ve£h)\\T7 


^r'8ln2  6Ii 


;  <ki  *kl> 

sin  0^  sin  6^ 


(4.66) 


For  convenience,  the  above  terms  are  summarized  below 


e  ■  Polarization  of  the  incident  electric  field  (E.  =  E*  e  ) 
a  i  o  a 

e.  =  Polarization  of  the  scattered  electric  field  (E  -  Es  e.  ) 
b  sob 


-*u*-h,y** 

/l  +  r  2  +  r  2~ 
A  Hx  Hy 


*►  A  A  A  A  A 

k^  *  k^^  ;  ^  -  -  sin  0.^  cos  ^  x  -  sin  0^,  sin  ^y-  cos  0.^  z 

k.'  «  k  k,1  ;  k.'  «  -  sin  0  cos  <b  x  -  sin  0  sin  <t>  y  -  cos  0  z 

ioi  i  s  s  s  s  s 


fk  A 

j  kiXn& 

"Ah  if>  w  a  i 


lci  *s*i  lh  l£lx"*l 


*;*"*  -i  -  - 

77-FT  eHv-ethXki 


Al *  A  A, 

*lv  *  ethxki 


.ln9tl-  Ikj-Sjl 


i£;xSni 


-*  A  A  ,  _  A 

r^-xx  +  yy  +  ^z 
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2e  cos  01 . 
r  £i 


Attempts  have  been  made  to  compare  (4.66)  with  the  equivalent  factor 
resulting  from  the  "tilted-plane"  approach  [5]  but,  unfortunately,  the  cor¬ 
respondence  is  not  easy  to  establish.  It  appears  that  such  a  comparison 
might  best  be  accomplished  by  comparing  numerical  values  of  (4.66)  with  the 
corresponding  factor  from  the  "tilted-plane"  approach  [5]. 

Equation  (4.65)  is  the  desired  result.  From  this  expression,  one  can 
easily  obtain  o  *  from 


lim  lim  f  4TtR2  1 
R-+00  A-*»  ^  A  E  2 
o 


(4.67) 


along  with  the  development  given  in  [6]  and  as  corrected  in  Section  2.  The 
contribution  of  the  zeroth  order  incoherent  power  <|6° Eab|  >  has  been  previ 
oualy  obtained  by  Sancer  [11]  and  his  results  can  be  used  directly  in  (4.67). 
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5.0  A  USEFUL  RELATIONSHIP  FOR  THE  JOINT  SLOPE  PROBABILITY  DENSITY  FUNCTION 


5.1  Background 

The  incoherent  power  scattered  in  and  about  the  specular  direction  de¬ 
pends  upon  the  joint  probability  density  function  (jpdf)  of  the  large  scale 
slopes.  The  jpdf  for  the  large  scale  slopes  is  also  important  in  determin¬ 
ing  the  degree  of  "tilt"  or  k-space  broadening  imparted  to  the  small  scale 
Bragg  scatterers.  Ideally,  one  would  like  to  measure  the  jpdf  for  the  large 
scale  slopes,  the  roughness  spectrum  of  the  small  scale  heights,  and  the  com¬ 
plex  dielectric  constant  of  the  surface  in  order  to  predict  the  average  scat¬ 
tering  properties  of  a  specified  section  of  terrain.  That  is,  these  surface 
measurements  would  be  substituted  in  the  rough  surface  scattering  model  which, 
in  turn,  would  provide  an  estimate  of  the  average  coherent  and  incoherent 
scattered  power.  From  a  practical  point  of  view,  measurements  of  the  jpdf 
for  the  slopes  and  the  small  scale  roughness  spectrum  are  very  difficult  to 
obtain  and  the  difficulty  increases  as  the  radar  or  electromagnetic  wave¬ 
length  decreases.  For  example,  in  the  case  of  an  L-band  system  with 

*  30  cm,  the  small  scale  part  of  the  scattering  model  will  require  sur¬ 
face  height  spectral  measurements  of  surface  undulations  having  wavelengths 
of  less  than  about  90  cm  because  *  (X  /2)csc  6.  for  backscatter. 

For  the  large  scale  features  of  the  surface,  the  jpdf  for  the  slopes  repre¬ 
senting  surface  features  having  spatial  wavelengths  greater  than  about  90  cm 
Is  required.  Obviously,  spectral  information  on  the  small  scale  features  is 
going  to  be  the  most  difficult  to  obtain.  However,  even  the  jpdf  for  the 
large  scale  slopes  is  going  to  be  difficult  to  estimate.  It  is  not  unrea¬ 
sonable  to  expect  that  we  can  obtain  measurements  of  the  jpdf  for  the  large 
scale  heights  and  even  the  correlation  function  for  the  heights  at  least  some¬ 
what  close  to  the  90  cm  spatial  resolution.  However,  this  information  must 


* 


somehow  be  translated  into  the  jpdf  for  the  surface  slopes  and  this  is  where 
the  difficulty  comes  in.  We  ignore  the  possibility  of  a  direct  measurement 
of  the  jpdf  of  the  large  scale  slopes  for  arbitrary  terrain  because  such  a 
task  appears  to  be  too  difficult  to  even  contemplate. 

The  question  basically  boils  down  to  the  feasibility  of  translating  or 
converting  measurements  of  the  jpdf  for  the  surface  heights  into  the  jpdf  for 
the  surface  slopes.  The  purpose  of  this  section  is  to  point  out  an  analyti¬ 
cal  means  for  accomplishing  this  transformation  and  suggest  that  the  scheme 
be  attempted  on  an  experimental  basis.  The  relationship  is  not  new  and,  in 
fact,  results  from  some  earlier  rough  surface  scattering  analysis.  However, 
it  has  apparently  gone  unnoticed  at  least  insofar  as  it  applies  to  this  very 
real  world  problem  of  translating  the  height  jpdf  into  the  slope  jpdf. 

5.2  The  Transformation 

Perhaps  the  oldest  approach  to  estimating  the  quasi-specular  incoherent 
power  scattered  by  a  rough  surface  is  now  called  the  autocorrelation  approach. 
Basically,  one  assumes  the  validity  of  physical  optics,  interchanges  the 
order  of  spatial  integration  and  ensemble  averaging  in  the  expression  for  the 
scattered  power,  and  assumes  Gaussian  surface  statistics  with  the  final  re¬ 
sult  that  the  average  scattered  power  is  dependent  upon  the  behavior  of  the 
surface  height  correlation  function  near  | Ar |  =  0[1]  .  In  the  mid-60's, 

Kodis  [2]  showed  that  the  average  scattered  power  could  alternatively  be 
interpreted  in  terms  of  the  number  of  specular  points  on  the  surface  and  the 
absolute  radii  of  curvature  at  the  specular  points.  Barrick  [3]  subsequently 
linked  these  two  approaches  in  the  high  frequency  limit  where  both  are  valid. 

In  the  process  of  establishing  the  similarity  between  the  authcorrela- 
tion  and  specular  point  approaches,  Barrick  obtained  a  relationship  between 
the  jpdf’s  of  the  surface  heights  and  the  surface  slopes.  In  particular. 


i 


66 


if  Pr  r  (C„»5  )  is  the  jpdf  of  the  x  and  y  surface  slope  components  and 

‘’x'.y  x  y 

4_  _  (k  ,k  )  is  the  joint  characteristic  function  for  the  surface  heights 
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(5.1) 


where  Ax  =  and  Ay  =  y^  -  y^  .  In  (5.1)  the  quantities  q^.q^  and 

q  are  limited  as  follows;  jq  |  £  1  ,  i  =  x,y,z  .  It  should  be  noted  that 
both  4  (•)  and  P  (•)  are  implicit  functions  of  the  surface  height 

ScS 

correlation  function;  this  is  how  the  (Ax, Ay)  variation  comes  about  in  (5.1). 

If  an  analytical  form  for  4_  is  available  then  (5.1)  can  be  used 

?1?2 

directly  to  obtain  the  jpdf  for  the  slopes.  In  cases  where  4_  (•)  is 

obtained  from  measured  data,  it  is  not  immediately  obvious  that  (5.1)  is  of 
any  practical  use  since  the  behavior  of  the  joint  height  characteristic  func¬ 


tion  will  not  be  known  in  the  limit  of  k  and  k  -m»  .  However,  consider 

x  y 

the  following  reasoning  as  a  means  for  obtaining  estimates  of  P  (*»•)  . 

^x^y 


Since  4  is  the  two-dimensional  Fourier  transform  of  the  jpdf  for 

the  height,  it  can  be  obtained  numerically  by  using  a  Fast  Fourier  Transform 
(FFT)  on  the  measured  jpdf  height  data.  The  result  of  this  operation  will  be 

A 

denoted  by  i,  .  Because  of  measurement  noise  and  particularly  quantiza- 

51>2 

A 

tion  noise  in  the  measured  height  Jpdf  data,  4  will  be  limited  to  values 

less  than,  say,  k  <  K  and  k  <  K  .  The  maximum  value  of  k  that  can 
’  x  ~  x  y  —  y  o 

be  achieved  in  (5.1)  is  therefore  max(K  /q  ,K  /q  )  .  If  q  is  small  then 

x  z  y  z  z 

the  resulting  maximum  value  of  kQ  can  be  very  large.  The  transform  variables 
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in  (5.1)  will  be  given  by  K  q  /q  and  K  q  /q  which  may  also  be  large, 

x  x  z  y  y  z 


depending  upon  q  and  q 
x  y 

puted  using  the  following; 


Thus,  if  q  is  near  zero  and  ?-  -  is  com- 
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it  may  turn  out  that  P  is  a  sufficiently  good  estimate  of  P  as  to 

be  useful  in  the  scattering  model.  Unfortunately,  this  approach  breaks  down 


when 


q^  or  q  •  0  ;  however,  it  may  be  possible  to  get  close  enough  to 


q^  *  0  or  q  ”  0  to  infer  the  behavior  of  V  along  these  lines  in  the 

X  y  ^xS 

qx»qy  -plane.  The  limits  on  the  integrations  in  (5.2)  symbolically  denoted 
as  ±  L  ,  will  be  determined  by  the  correlation  length  of  the  surface,  i.e. 
the  separation  distance  for  which  the  surface  height  correlation  function 
is  essentially  zero. 

An  alternate  approach  to  estimating  P  is  to  examine  the  asymptotic 

^x^y 

behavior  of  _  as  k  -+K  and  k  -+-K  .  From  this  behavior,  it  may  be 
SjC2  x  x  y  y 

possible  to  generate  an  asymptotic  functional  dependence  of  9  on  k 

hh  * 


and  k 


By  repeating  this  procedure  for  different  values  of  Ax  and  Ay  , 


it  might  be  possible  to  also  generate  or  build-in  the  functional  dependence 

A  A. 

of  9  on  Ax  and  Ay  .  In  this  manner,  the  dependence  of  $  upon 

^1^2  Cl4 

k  ,k  ,A  ,  and  A  is  obtained  at  least  in  the  limit  of  moderately  large  k 
x  y  x  y  x 

and  ky  .  This  functional  form  could  then  be  transformed  according  to  (5.2). 
The  major  problem  here  is  that  the  accuracy  of  the  result  will  depend  directly 
upon  how  precisely  the  surface  height  correlation  function  is  known  near 
Ax*0  and  Ay » 0  .  This  statement  results  from  the  fact  that  the  behavior 


68 


of  the  transform  of  a  function  as  k-»°°  is  directly  determined  by  the  behavior 
of  the  function  as  Ar  +  0  [4]. 

The  problem  of  converting  height  jpdf  data  into  slope  jpdf  results  is 
definitely  not  easy.  Even  with  the  use  of  (5.1)  the  problem  still  poses  a 
number  of  numerical  complexities,  primarily  because  of  the  required  limit  as 
kQ  +00  .  However,  as  discussed  above,  (5.1)  does  provide  some  hope  in  solving 
what  is  otherwise  a  totally  untractable  problem.  It  is  felt  there  is  suf¬ 
ficient  hope  as  to  warrant  further  investigation  of  the  utility  of  (5.1)  In 
the  solution  of  this  problem. 
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